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ABSTRACT 


A description  of. a linear  model  for  forecasting  short 
term  utility  loads  is  presented.  The  model  is  used  to 
simulate  the  improvement  to  the  utility  load  factor  by 
providing  a basis  for  scheduling  deferred  loads. 

Simulations  were  performed  using  data  from  Central  Vermont 
Public  Service  Corporation  for  the  period  1976-1977.  The 
results  indicate  that  an  average  hourly  prediction  error  of 
less  than  5%  is  possible  with  this  model,  and  that  an 
increase  in  the  utility  system  load  factor  is  directly 
related  to  the  magnitude  of  the  controllable  load. 
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CHAPTER  I 


INTRODUCTION 


I.  Introduction 

An  electric  utilities'  loads  are  graphically  examined  and  a 
basic  terminology  is  presented. 

I.  1.0  Objectives 

The  primary  purpose  of  this  study  is  to  investigate  the 
impact  which  scheduling  of  deferred  loads  will  have  on  the  load 
factor  of  utilities.  Collaterally  the  problem  of  short  term  load 
forecasting,  in  order  to  optionally  schedule  loads,  became  an 
equally  important  area  for  investigation.  Once  the  future  load 
pattern  has  been  forecast,  then  for  a given  controllable  load  the 
off  peak  system  load  be  increased  with  potential  benefits  to  bxjth 
the  utility  and  the  consumer.  These  features  were  siniulated 
using  a two  year  period  of  load  data  from  a Vermont  utility. 

I.  2.0  Background 


A basic  understanding  of  the  electric  utility  industry  in  the 
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United  States  and  particularly  in  New  England  was  required  as  a 
starting  point  for  the  load  forecasting  phase.  Prior  knowledge 
of  system  loads,  a load  forecast,  would  enable  a utility  to 
schedule  its  generating  capacity  to  match  its  loads  and  hence 
reduce  generating  costs.  In  actuality  individual  utilities  in 
Vermont  have  no  such  short  term  prediction  available.  Due 
primarily  to  the  emergence  of  power-pools  in  the  mid  1960 's 
utilities  no  longer  scheduled  the  operation  of  their  own 
generating  systOTS.  This  function  was  assumed  in  a regional 
basis  by  the  power  pools  who  in  turn  developed  forecasting 
procedures  for  a much  larger  geographic  area  and  higher  system 
load.  The  number  of  technical  papers  published  on  load 
forecasting  seens  to  have  peaked  in  the  early  1970 's  (Sachdev, 
1977) . 

I.  3.0  Utility  Responsibilities  and  Definitions 

/ 

Most  electric  utilities  are  private  corporations  regulated  by 
local  governmaits  in  the  public  interest.  Until  recently  each 
has  been  a monopoly  within  its  area  of  responsibility.  By  law 
the  utility  is  required  to  reliably  supply  all  of  the  power  its 
customers  demand,  where  power  is  defined  as  the  electric  load,  or 
total  utility  system  load,  at  time  t.  It  is  generally  integrated 

? 

0 

I 


3 


over  an  hour  and  recorded  in  units  of  megawatt  hours  (MWH) . The 
demand  varies  over  each  time  interval  examined  from  seconds 
through  years.  Short  term  load  fluctuation  defined  as  those 
occurring  from  one  through  twenty- four  hours  in  the  future  are 
important  to  the  system  operator,  the  individual  responsible  for 
insuring  the  uninterrupted  flow  of  power.  Longer  term 
forecasting  does  not  imnediately  apply  to  the  system  operator  and 
often  involves  socio-economic  factors.  Shorter  term  forecasting 
(i.e.  those  system  load  fluctuations  occurring  fran  seconds  to 
several  minutes)  appear  to  be  white  noise  and  are  almost  always 
small  enough  to  be  balanced  by  automatically  adjusting  the 
operating  frequency  of  generating  units  (Galiana,  1971) . It  has 
been  estimated  that  a 12%  imbalance  between  generation  and  load 
results  in  a frequency  change  of  1 Hz.  Major  blackouts,  such  as 
those  experienced  by  New  York  City,  are  characterized  by  line 
frequency  shifts  of  from  one  to  two  Hz  (Ewart,  1978) . 

/ 

The  term  'load  management'  is  defined  as  any  action  taken  by 
an  electric  utility  to  manage  its  loads.  Load  managenent  is 
subdivided  into  supply  management  where  utilities  modify  the 
system  load  as  seen  from  the  generator  ports  without  affecting 
individual  customer's  consumption  patterns  and  customer 
management  where  the  power  demand  pattern  itself  is  modified  in 
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sane  manner.  The  latter  load  management  strategy  cocimonly  takes 
the  form  of  promotion  of  conservation  measures,  promotion  of  off 
peak  electric  periods,  or  a rate  structure  v^ich  penalizes  peak 
period  consimption.  Only  recently  in  the  United  States  has 
direct  load  control  by  the  utility  been  added  to  the  available 
methods  of  customer  load  management  (Galiana,  1977) . Here  the 
utility  actually  controls  or  schedules  the  pattern  of  consumption 
in  some  manner.  The  operation  of  selected  appliances  may  be 
controlled  by  radio  frequency  signals,  coded  pulses  superimposed 
on  the  60  Hz.  line  voltage  (i.e.  ripple  control) , or  load  timers. 
For  the  first  time  the  system  operator  can  shape  the  daily  load 
curve  and  selectively  shed  loads  in  case  of  system  anergencies. 


Shaping  the  load  curve  to  improve  the  system  load  factor 
impacts  directly  on  the  cost  of  power  to  the  consumer.  The  load 
factor  is  a measure  of  the  efficiency  of  the  electrical  system. 
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It  is  defined  by  the  industry  as: 


j^max  (t)j  X 24 

Z = Systan  Load  MWH 
n = Time  Periods 
t = Time 

Normally  a utility  measures  its  efficiency  by  its  yearly  load 
factor  i.e.  n = 8760  hours/year.  Central  Vermont  Public  Service 
Utility's  1977  load  factor  was  approximately  60%.  Obviously  the 

simplest  way  to  improve  load  factor  is  to  reduce  the  maximum 

/ 

systan  demand.  Controlling  a substantial  load  through  ripple 
signals  and  then  being  able  to  defer  or  shift  that  load  to  an  off 
peak  period  accomplishes  this  reduction  in  maximum  ^ (t  ) . 

Deferrability  of  most  electric  loads  from  the  consumer's 
viewpoint  is  undesirable.  The  customer  wants  his  lights  to  be  on 
and  appliances  to  work  when  he  closes  the  switch.  Similarly 
industrial  motor  loading  is  non-deferrable.  However  substantial 
electric  power  is  consumed  in  heating  homes  and  water.  For  a 
Vermont  electric  home  the  average  heating  load  was  19,000  KWH 
from  a total  annual  load  of  25,000  KWH  (Vermont  Load  Study, 
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1974) . Resistive  heating  loads  are  deferrable  to  a certain 
extent  since  this  electric  energy  can  be  stored  in  a heated 
medium  for  later  consumption.  Thermal  Energy  Storage  (TES)  is 
the  commonly  used  term  for  this  process.  In  Germany,  buildings 
have  employed  TES  heating  for  over  twenty  years  and  in  certain 
regions  it  now  accounts  for  over  45%  of  the  region's  peak  load. 
There  are  now  more  than  six  million  electrically  heated 
residences  throughout  the  United  States,  so  the  potential  for 
load  factor  improvement  through  TES  and  customer  load  managanent 
exists.  It  is  not  a completely  new  or  untried  technology 
(Siemens,  1978) . 

Energy  costs  in  Vermont  are  more  than  that  of  any  other  state 
in  the  continental  United  States.  The  weighted  average  price  of 
one  million  fossil  fuel  BTUs  paid  by  utilities  was  61.7  cents  in 
1973.  Today  that  same  figure  is  $2.30  and  in  1973  dollars  it  is 
still  $1.80.  That  is  a 373%  increase  in  four  years  and  also 
means  that  Vermont  utilities  pay  780%  as  much  for  fossil  energy 
as  the  nation's  cheapest  state  - North  Dakota  (Buggraf,  1978). 
Making  the  energy  system  more  efficient  is  one  way  of  reducing 
this  cost  disparity. 
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I.  4.0  Load  Patterns 

The  question  of  forecasting  what  the  system  load  will  be  has 
been  approached  in  the  literature  many  ways;  however,  all 
forecasting  has  been  in  some  measure  based  on  the  total 
historical  load  pattern  as  opposed  to  a decor^xjsition  of  the 
load.  Hence  the  following  pattern  analysis  will  be  a starting 
point.  Data  from  Central  Vermont  Public  Service  Corp.  was  chosen 
due  to  its  availability,  although  data  from  any  utility  could 
have  been  used. 

I.  4.1  Monthly  Pattern 

Fig.  1 illustrates  the  periodic  behavior  of  Central  Vermont 
Public  Utility's  load  for  the  month  of  April  1976.  Note  the 
weekend  valleys  and  the  general  inverse  load  trend  with 
increasing  average  temperature  throughout  the  period.  This 
perspective  of  the  load  was  chosen  in  oraer  to  prevent  the 
morning  peak  from  masking  the  noon  valley  and  evening  peak. 
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I.  4.2  Weekly  Load  Patterns 

Fig.  2 illustrates  the  system  load  for  a Sunday  through 
Saturday  period  in  May  1976.  Again  note  the  similarity  of 
weekday  versus  weekend  load  levels.  The  pattern  is  especially 
pronounced  on  Tuesday  through  Friday  and  has  lead  some  authors  to 
consider  the  normal  system  load  problem  as  being  made  up  of  three 
regions  i.e.  weekends,  Tuesdays  through  Fridays  and  Mondays 
(Wernhof,  1978). 

I.  4.3  Repeated  Weekday  Pattern 

A deeper  graphical  appreciation  of  the  similarities  in  system 
load  for  weekdays  can  be  gained  from  Figures  3 and  4 vdiich  show 
load  on  consecutive  Mondays  and  Thursdays  in  May  1977.  Mondays 
lack  the  pronounced  noonday  load  valley.  Note  the  affect  which 
the  Memorial  Day  holiday  had  on  the  load.  Holidays,  of  course, 
must  be  forecasted  on  a special  case  basis.  A pronounced  noon 
valley  is  characteristic  of  many  utility  systems  and  gives  rise 
to  the  phrase  "camel's  hump  load  pattern".  Note  eilso  that  the 
load  does  not  fluctuate  greatly  from  one  week  to  the  next.  Even 
in  the  seasonal  translation  months  the  total  daily  load  did  not 
change  by  more  than  20%  over  a weeks  time. 


9-15  1976  Loads 
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I.  4.4  Similar  Days  Separated  by  N Periods 

Fig.  5 shows  the  system  loads  for  the  second  week  of 
September  1976  and  the  same  week  one  year  later.  The  basic 
pattern  is  of  course,  the  same;  however,  the  electrical  load  has 
evolved  enough  to  preclude  merely  using  the  "same  day  last  year" 
as  a basis  for  prediction  this  year.  A legitimate  question  here 
would  be,  "Is  a similar  day  last  year  a good  prediction  for  this 
year?".  Neglecting  until  Section  III  a discussion  of  what 
constitutes  a good  prediction  leads  to  the  problem  of  defining  a 
similar  day.  Similar  in  v^at  manner?  Date,  day  of  week, 
maximum  and  minimum  temperature  range,  wind,  cloud  cover, 
precipitation,  and  weather  trends  are  all  possible  factors  in 
characterizing  similarity.  There  is  a high  probability  that 
similarity  exists  only  in  a subjective  manner,  especially  since 
the  system  peak  day  is  unique  each  year,  and  it  has  occurred 
randomly  in  the  sixty  day  December-January  time  window.  Pig.  6 
illustrates  the  similarity  of  Thursdays  separated  by  one  week, 
one  month,  and  one  year.  Their  average  daily  temperature  was 
within  5 degrees  Fahrenheit  throughout  to  minimize  the 
temperature  dependence  of  load.  Temperature  data  for  the 
utilities  visited  was  maintained  in  Fahrenheit,  and  there  was  no 
discussion  of  plans  to  convert  to  Centegrade. 


Fig.  5 Loads  for  Week  of  Sept.  1976  & 1977 
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I.  4.5  Peak  Day  Pattern 

Fig.  7 shove  the  load  pattern  for  the  Central  Vermont  Public 
Utility  system  peak  day,  Tuesday,  18  January  1977.  Also  shovm  is 
Tuesday  16  August  1977  to  illustrate  the  seasonal  load  variation 
in  total  daily  load  frcm  8036  in  January  to  4987  MW  in  August. 
Average  temperature  on  the  peak  day  was  -5*  F versus  69*  F on  16 
August  1977.  Vermont  is  definitely  a winter  load  peaking  state 
as  opposed  to  the  New  York  megalopolis  which  has  an  air 
conditioning  suitmer  peak. 

The  variation  in  load  from  season  to  season  is  in  general  due 
to  the  length  of  day  and  weather,  and  this  factor  in  turn 
directly  impacts  on  the  customer's  consumption  patterns.  For 
example,  in  the  winter  electric  heating  and  indoor  cooking 
increase  as  do  lighting  loads  due  to  shorter  days.  Plant 
shutdowns  for  vacation  are  rare  in  the  winter  and  combined  with 
their  need  for  electric  heating  increase  the  industrial  system 
load.  Alternatively  during  the  sumter  refrigeration  loads 
increase,  lighting  loads  decrease,  and  people  generally  spend 
less  time  inside.  Industrial  loads  decrease  due  to  vacation 
schedules.  Sunmer  peaking  will  not  occur  in  Vermont  because  the 
ski  industry  and  heating  consume  so  much  electricity. 


Low  System 
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CHAPTER  II 

ANALYSIS 

II.  1.0  Introduction 

The  method  of  load  data  collection  and  other  available  data 
sources  are  discussed.  Standard  statistical  analysis  techniques 
are  applied  to  system  load  data  to  determine  which  variables  are 
significant.  System  load  relationships  with  temperature,  (i.e. 
lag  and  proportionality)  are  examined. 

II.  1.1  Utility  Load  Data 

Data  for  this  study  was  made  available  by  Central  Vermont 
Public  Service  Corporation  at  their  Rutland  headquarters.  The 
systan  dispatcher  maintains  a handwritten  hour  by  hour  log  of 
system  loads,  weather  in  Rutland,  and  operating  status  on  desk 
sized  ledger  sheets.  These  in  turn  are  bound  and  constitute  the 
legal  historical  operating  record  of  the  utility.  Two  years  of 
this  data  was  copied,  formatted  and  entered  on  the  University  of 
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Vennont  Canputer  System.  The  Julian  Data  Numbering  System  was 
chosen  to  organize  the  individual  daily  files  of  hourly  data 
(Appendix  1).  The  entire  data  set  consisted  of  over  17,500 
records  organized  by  day  as  shown  (Appendix  1) . This  data  had 
not  previously  been  available  for  ccnputer  analysis. 

II.  1.2  Weather  Forecasts 

While  not  all  load  prediction  schones  in  the  literature  use 
weather  forecasts,  all  actual  forecasting  employs  a weather 
factor.  The  accuracy  of  weather  forecasts  will  be  a limiting 
factor  in  any  load  prediction  scheme.  The  National  Weather 
Service  is  able  to  provide  twenty-four  hour  forecasts  of  the 
maximum  and  minimum  tatiperatures  within  two  degrees  Fahrenheit 
more  than  95%  of  the  time.  Their  prediction  will  include  a 
temperature  profile  which  indicates  vrtiether  there  will  be  a 
morning  low  and  afternoon  high  or  the  converse  shape. 
Approximately  90-95%  of  the  days  of  a year  follow  the  morning  low 
evening  high  temperature  pattern.  Specific  hourly  temperature 
forecasts  are  not  currently  available.  While  commercial  weather 
forecasting  bureaus  claim  to  provide  a more  accurate  or  'better' 
forecast  the  experience  of  New  England  Power  Exchange  (NEPEX)  has 
not  supported  this  claim.  NEPEX  no  longer  subscribes  to  a 


20 


private  weather  service.  Other  parameters  such  as  wind  or  cloud 
cover  are  not  available  on  an  hourly  basis  with  any  degree  of 
accuracy.  A major  limiting  factor  in  load  forecasting  is  the 
quality  and  quantity  of  weather  information  since  all  that  is 
reliably  available  is  a forecast  of  tomorrow's  maximum  and 
minimum  temperature. 

II.  1.3  Consumption  Patterns 

Extensive  work  has  been  done  (The  Vermont  Load  Study,  1975) 
catagorizing  consunption  patterns  of  various  components  of  the 
electrical  system  load.  Germane  to  this  study  is  the  result  that 
power  consumption  in  electrically  heated  homes  slightly  leads  the 
total  system  load  throughout  a twenty- four  hour  period.  Fig.  8 
illustrates  this  point  and  also  shows  that  the  heating  load  has 
no  dead  spots  or  zero  hours.  On  the  average  this  allows  the 
system  load  curve  to  be  decreased  by  some  constant  factor  for 
each  hour  of  controlled  operation. 

A 1975  study  (Veno,  1975)  of  the  economic  viability  of  TES 
units  in  the  Central  Vermont  Public  Utility  areas  concluded  that 
a peak  day  load  factor  improvement  approaching  18%  was  possible 
by  scheduling  TES  units  from  8pm  to  6am. 


X 
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VERMONT  LOAD  STUDY 
GROUP  II  - RESIDENTIAL 
(SPACE  HEATING  CUSTOMERS) 

AVERAGE  WEEKDAY  HOURLY  LOADS  - JANUARY  1975 


Fig.  8 
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In  that  study  no  provision  was  made  to  forecast  the  loads  or 
consider  varying  hourly  heating  consumption.  A constant  MW  value 
determined  by  the  number  of  electric  heating  customers  and 
average  installed  heating  capacity  was  subtracted  from  the 
observed  peak  day  hourly  load.  The  corresponding  TES  load  was 
then  scheduled  during  the  nighttime  valley.  Fig.  9 illustrates 
the  procedure.  The  concept  of  subtracting  an  average  load  due  to 
electric  heating  from  the  daily  load  curve  will  also  be  anployed 
in  this  study. 

II.  2.0  Data  Analysis 

A series  of  analyses  were  performed  on  data  files  using  the 
Statistical  Package  for  the  Social  Sciences  (SPSS) . The  first 
objective  was  to  determine  the  relative  strength  of  relationships 
among  the  variables,  and  the  second  was  to  employ  standard  linear 
regression  analysis  as  a forecasting  technique. 


TES  Procedure 
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II.  2.1  Correlation 

Initial  statistical  tests  were  performed  for  a given  month  of 
the  year,  day  of  the  week  and  hour  of  the  day,  considering  the 
following  as  potential  variables; 

System  Load 

Tanperature  Fahrenheit 
Cloud  Cover 
Wind  Speed 

For  the  month  of  October  1977  scattergrams  were  obtained  relating 
system  load  to  temperature  tor  all  30  days.  There  was  no 
discernible  trend.  For  the  same  data  a scattergram  relating  MW 
to  hour  of  the  day  produced  a band  whose  shape  resanbled  a daily 
load  curve. 
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The  Kendall  and  Spearman  Correlation  Coefficients  were  generated 
for  October  1977  system  load  (MW)  with  tonperature,  wind,  cloud 
cover,  and  hour  of  the  day  as  variables: 


M^-TEMP 

MW-HR 

MW-CLOJD 

MW-WIND 

Kendall 

.1378 

.2686 

-.0041 

J582 

Spearman 

.2176 

.4567 

-.0089 

.2088 

Significance 

.001 

.001 

0.4 

.001 

TABLE  1 Correlation  Coefficients 

Both  correlation  methods  show  that  load  is  most  dependent  on  hour 
of  the  day,  however  vAiether  wind  or  temperature  is  the  second 
strongest  independent  variable  is  not  clear.  The  correlation  of 
load  with  cloud  is  not  significant  probably  due  to  the  weather 
coding  scheme  employed  for  cloud  cover  which  reduces  a large 
range  of  data  into  a relatively  small  number  of  catagories 
(Appendix  1) . 
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The  Pearson  Correlation  Coefficients  for  October  1976,  are 
presented  in  Table  2. 


MW 

HOUR 

TEMP 

CLOUD 

1.0 

.4823 

,0344 

.0254 

HOUR 

.4823 

1.0 

.175 

-.034 

TEMP 

.0344 

.175 

1.0 

.054 

CLOUD 

.0254 

-.0342 

.054 

1.0 

WIND 

.106 

.0159 

.008 

.1123 

TABLE  2 Pearson  Correlations 


II.  2.2  Hourly  Correlations 

One  would  specifically  expect  an  inverse  correlation  between 
temperature  and  system  load  during  the  winter  in  a state  like 
Vermont.  To  produce  a valid  relationship  between  system  load  and 
temperature  it  is  necessary  to  remove  life  cycle  affects  from  the 
data.  Hence  each  of  the  twenty- four  hours  was  treated  as  a 
separate  variable  and  correlations  performed  for  that  hour 
between  MW  and  temperature.  Results  of  hourly  correlations  for 
each  of  six  monthly  time  periods  are  presented  in  Table  3. 
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2 

Correlation  squared  (r  ) indicates  the  proportion  of 

variation  in  system  load  explained  by  temperature  (i.e.  for 

2 

October  1976  at  lAM  r = -.618,  r = .382) , and  thus  38%  of  the 
variation  in  load  is  explained  by  linear  regression  on  the 
tanperature  variable  alone.  All  coefficients  listed  are 
significant  at  the  .001  level. 

Examining  partial  correlations,  or  the  relationship  between 
two  variables  vrfiile  adjusting  for  the  effects  of  other  variables, 
led  to  the  conclusion  that  for  October  1977  system  loads,  the 
hour  of  the  day  is  the  predominant  independent  variable  followed 
by  tanperature,  then  much  weaker  cloud  cover  and  finally  wind 
speed.  In  the  regression  analysis  performed  the  r ^ change  due  to 
cloud  was  less  than  .002  vrtiile  that  for  wind  was  less  than  .0005. 
The  load  forecaster  at  NEPEX  had  stated  that  cloud  cover  was  more 
significant  than  wind  to  his  system  load  (Wernhof , 1977) . 
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REGRESSION 


OCTOBER  1976  OCTOBER  1977 


HOUR 

COEFFICIENTS 

COEFFICIENTS 

A 

B 

CORRELATION 

A 

B 

CORRELATION 

1 

197.15 

-.982 

-.618 

135.6 

-.427 

-.255 

2' 

184.31 

-.896 

-.591 

166.78 

-.235 

177 

3 

181.08 

-.916 

-.633 

174.5 

-.503 

-.503 

4 

182.49 

-.970 

-.632 

170.3 

-.493 

-.401 

5 

189.0 

-1 . 125 

-.699 

190.92 

-.918 

-.639 

6 

207.64 

-1.26 

-.79 

193.22 

-.624 

-.568 

7 

249.35 

-1.364 

-.821 

240.5 

-.842 

-.603 

8 

314.31 

-1.59 

-.889 

309.4 

-1.08 

-.584 

9 

326.3 

-1.37 

-.898 

330.  18 

-1.05 

-.754 

10 

327.24 

-1.33 

-.859 

314.04 

-.739 

-.389 

11 

318.23 

-1.135 

-.932 

358.35 

-1.63 

-.782 

12 

319.23 

-1.21 

-.857 

358.95 

-1.63 

-.638 

13 

302.99 

-1.08 

-.839 

345.6 

-1.57 

-.805 

14 

305. 12 

-1.00 

-.814 

336.84 

-1.37 

-.712 

15 

300.83 

-1 . 107 

-.85 

327. 15 

-1.303 

-.830 

16 

286. 17 

-.973 

-.785 

328.04 

-1.396 

-.858 

17 

293.35 

-1 .084 

-.784 

325.  16 

-1.29 

-.787 

18 

307.76 

-1.114 

-.761 

343.18 

-1.51 

-.743 

19 

313.84 

-1.225 

-.844 

319.56 

-1 .05 

-.634 

20 

305.0 

-1  .08 

-.777 

283.27 

-.545 

-.333 

21 

290.84 

-.9939 

-.753 

269.76 

.372 

-.165 

22 

276.79 

-1.02 

-.742 

286.82 

-.958 

-.508 

23 

257.24 

-1.09 

-.749 

254.52 

-.797 

-.530 

24 

224.87 

-1.02 

-.805 

213.8 

-.458 

-.299 

TABLE  3 Regression  & Correlation 
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JANUARY  1976  FEBRUARY  1976 


HOUR 

COEFFICIENTS 

COEFFICIENTS 

A 

B 

CORRELATION 

A 

B 

CORRELATION 

1 

260.55 

-1.68 

-.796 

242.55 

-1.45 

-.34 

2 

248.05 

-1.7 

-.85 

233.9 

-1.54 

-.38 

3 

242. 13 

-1.69 

-.865 

234.5 

-1.7 

-.90 

U 

240.66 

-1.61 

-.834 

227.33 

-1.55 

-.89 

5 

238.9 

-1.71 

-.84 

228.74 

-1.53 

-.896 

6 

250.66 

-1.76 

-.79 

217.6 

-.566 

-.41 

7 

279.3 

-1.5 

-.57 

263.0 

-1.29 

-.53 

3 

313.0 

-.77 

-.31 

284.96 

-.558 

-.229 

9 

325.7 

-.55 

-.25 

327.5 

-1.48 

-.496 

10 

337.1 

-1.03 

-.41 

336.97 

-1.61 

-.639 

11 

335.6 

-.87 

-.34 

309.7 

-.578 

-.28 

12 

332.99 

-.75 

-.30 

332.5 

-1.44 

-.65 

13 

318.32 

-.59 

-.28 

316.5 

-1.17 

-.67 

14 

313.8 

-.341 

-.  147 

310.78 

-1 . 15 

-.52 

15 

297.66 

-.017 

-.007 

303.4 

-1.14 

-.52 

16 

300.43 

-.372 

-.  175 

293.8 

-1 . 19 

-.53 

17 

316.76 

-.316 

-.  148 

308.24 

-1.23 

-.57 

18 

345.86 

-.765 

-.382 

339.4 

-1.78 

-.72 

19 

340.76 

-.836 

-.498 

331.'* 

-1.24 

-.53 

20 

339.4 

-1.11 

-.582 

331.44 

-1.49 

-.58 

21 

331.47 

-1 .24 

-.668 

323.75 

-1.63 

-.73 

22 

312.76 

-1.00 

-.598 

313.5 

-1.85 

-.82 

23 

307.57 

-1.84 

-.74 

291.9 

-1.74 

-.81 

24 

275.8 

-1.33 

-.706 

265.5 

-1.61 

-.83 

JANUARY 

1977 

FEBRUARY  1977 

HOUR 

COEFFICIENTS 

COEFFICIENTS 

A 

B 

CORRELATION 

A 

B 

CORRELATION 

1 

245.85 

-1.49 

-.614 

253.3 

-2.15 

-.878 

2 

234.6 

-.939 

-.46 

237.1 

-1.9 

-.87 

3 

228.8 

-1.19 

-.579 

231.7 

-1.87 

-.89 

4 

224.36 

-.56 

-.37 

230. 1 

-1.75 

-.77 

5 

225.3 

-.624 

-.479 

226.4 

-1.49 

-.72 

6 

233.6 

-.474 

-.418 

224.8 

-.65 

-.40 

7 

269.8 

-1.18 

-.40 

257.0 

-.728 

-.344 

8 

307.3 

-.21 

-.098 

295.9 

-.47 

-.183 

9 

327.2 

-.62 

-.175 

348.6 

-2.66 

-.66 

10 

323.0 

-.238 

-.069 

350.98 

-2.56 

-.62 

11 

339.2 

-1.59 

-.365 

356.5 

-2.54 

- . 606 

12 

334.4 

-1 .09 

-.285 

356.9 

-2.5 

-.653 

13 

324.95 

-1 . 16 

-.333 

341.6 

-2.2 

-.66 

14 

290.8 

f.752 

+ . 175 

356.6 

-2.8 

-.643 

15 

302.7 

-.389 

-.08 

337.7 

-2.38 

-.57 

16 

309.2 

-1 .07 

-.28 

331.6 

-2.29 

-.59 

17 

320.0 

-.81 

-.23 

341.4 

-2.5 

-.64 

18 

346.8 

-1.02 

-.285 

363.6 

-2.6 

-.72 

19 

341.6 

-.835 

-.278 

360.9 

-2.3 

-.72 

20 

332.7 

-.52 

-.31'< 

346.9 

-2.1 

-.72 

21 

326.04 

-1.03 

-.365 

335.6 

-2.02 

-.73 

22 

302.3 

-.196 

-.  133 

317.4 

-2.06 

-.795 

23 

288.2 

-.995 

-.38 

292.7 

-1.92 

-.727 

24 

262.6 

-1.12 

-.485 

264,4 

-1.7 

-.75 

TABLE  3 Regression  & Correlation 
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These  values  are  so  small  that  the  corresponding  variables  do 
not  affect  the  load.  Hence  as  an  initial  forecasting  methodology 
neither  cloud  cover  nor  wind  will  be  included  in  the  prediction 
routine . 

II.  2.3  Heating  Lag 

The  existence  of  a lag  between  temperature  change  and  system 
load  change  has  been  postulated  (Hajdu,  1977) . For  a temperature 
change  the  lag,  when  it  exists  is  between  four  to  eight  hours, 
and  is  attributed  to  heat  storage  or  inertia  in  massive  buildings 
such  as  would  be  common  in  large  city  loads. 

As  a method  of  identifying  the  lag  of  Central  Vermont  Public 
Utility's  load  five  Wednesdays  in  March  1977  were  plotted  in  Fig. 
10  and  Fig.  11.  While ‘the  inverse  relationship  between 
temperature  and  system  load  is  evident  there  does  not  appear  to 
be  a lag  affect.  Points  where  the  daily  temperature  profiles 
cross  closely  correspond  to  load  profile  crossings.  The 
conclusion  is  that  Centred  Vermont  does  not  have  a time  lag 
between  system  load  and  temperature  change. 
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Temperatures  March  Wednesdays  1977 
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II.  2.4  Temperature  Affects 

The  inverse  relationship  of  temperature  and  MW  can  be  clearly 
seen  in  Fig.  12  vrtiich  shows  both  the  hourly  systan  load  for  the 
second  week  of  March  1977  and  the  hourly  temperature.  By 
canparing  the  slope  of  the  temperature  with  the  increasing  load 
an  approximate  relationship  is  evident.  March  loads  increase 
approximately  one  megawatt  for  each  one  degree  drop  in 
temperature . 

Fuel  oil  companies  commonly  use  a degree  day  measure  to 
determine  v^en  additional  fuel  deliveries  are  required.  The 
degree  day  is  a measure  of  the  declination  of  the  mean  daily 
outdoor  temperature  from  an  arbitary  neutral  comfort  level  of  65 
degrees.  Positive  degree  days,  temperature  less  than  65*, 
generally  mean  that  energy  is  being  used  for  heating.  While  not 
commonly  used,  negative  degree  days  partially  measure  air 
conditioning  requirements,  neglecting  humidity.  Degree  days  are 
a post  hoc  measure  of  the  tonperature  profile,  and  because  of 
this  they  are  not  used  in  forecasting  electric  loads  (Wernhof, 


1978). 
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CHAPTER  III 


PREDICTION  METHODS 


III.  1.0  Introduction 

Several  methods  of  system  load  modeling  and  forecasting  are 
examined  and  these  methods  are  applied  to  develop  a model  for 
Central  Vermont  Public  Service  Corp.  A simple  hybrid  forecasting 
algorithm  is  developed  and  found  to  yield  acceptable  results,  and 
a computer  program  implanenting  the  algorithm  is  presented. 

III.  1.1  Catagories 

Historically,  when  individual  utilities  controlled  their 
generating  resources  and  made  their  own  daily  load  forecasts, 
the  prediction  was  based  on  years  of  operating  experience.  The 
method  then  used  to  arrive  at  a forecast  was  an  extrapolation  of 
past  system  load  experience  manually  correlated  or  modified  by 
the  weather  prediction.  The  correlation  criteria,  while 
subjective,  were  dynamic  and  could  be  complex.  Due  to  the  lack 
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of  a fixed  mathonatical  methodology,  consistent  forecasts  for  the 
same  criteria  were  not  possible.  This  did  not  mean  that  the 
operator  produced  forecasts  were  invalid.  In  fact.  Hydro  Quebec, 
which  is  reputed  to  have  one  of  the  best  prediction  schemes  in 
use,  relies  on  an  operator's  judgement  to  fine  tune  the  ccmputer 

1 

produced  forecast  for  weather  affects  {Bechard,  1978) . New 
England  Power  Exchange  (NEPEX) , a scheduling  entity  whose  peak 
hour  demand  is  approximately  16,000  MW,  relies  solely  on  the 
operator  technique  (Wernhof , 1978) . Results  obtained  with  an 
operator's  forecast  range  frcm  5 to  10%  average  error  over  a 
twenty-four  hour  period,  and  as  mentioned  earlier,  several 
utilities  which  had  experimented  with  mathematical  forecasts  have 
returned  to  the  operator  method  (Galiana,  1975) . 

The  forecasting  method  employed  should  logically  be  related 
to  the  type  of  forecast  required.  What  does  the  utility  need  to 
know  in  the  future?  Is  it  the  next  day's  hourly  load,  only  the 
peak  load,  total  daily  energy  supplied,  the  system  load  at  a 
specific  hour  for  system  maintenance  or  only  loads  for  security 
purposes?  The  same  criteria  which  answer  these  questions  should 
be  used  to  judge  the  results  of  a forecasting  method.  Judging  an 
hourly  forecasting  procedure  by  comparing  the  total  actual  daily 
load  with  the  sunmation  of  24  hourly  forecasts  is  not  an  adequate 
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indicator  of  hourly  accuracy  (Bechard,  1978) . 

Mathematical  nxx3els  for  the  load  are  divided  into  three  broad 
catagories  depending  on  the  inclusion  of  a weather  forecast 
variable  (Galiani,  1975).  These  are: 

1.  Models  vrtiich  include  weather  effects  but  do  not 

rely  on  the  latest  load  behavior  for  the  forecast. 

2.  Models  which  rely  on  time  of  day  and  latest  load  data. 

3.  Models  which  include  weather,  time  of  day  and  immediate 

past  load  in  the  forecast. 

Arguments  have  been  made  for  the  validity  of  each  model  type,  but 
by  inspection  Model  3 is  most  general  since  it  includes  the 
largest  nunber  of  available  data  inputs. 

III.  1.2  Time-Series  Models 

The  most  frequently  encountered  model  based  entirely  on  prior 
loads,  in  the  literature  on  load  forecasting,  is  the  time-series 
model.  Since  the  late  1960 's  the  Box-Jenkins  approach  has 
dominated  time-series  investigations.  A time-series  is  used  to 
fit  the  observations  of  system  load  made  at  discrete  hourly  time 
intervals.  The  method  relies  on  the  premise  that  the  load  curve 
has  a recognizable  shape  vrtiich  is  periodic  and  which  will  repeat 


38 


itself  in  the  future  to  fit  a function  of  time  to  the  load.  The 
time-series  approach  normally  does  not  afford  greater  weigh  to 
the  most  recent  load,  except  by  reestimating  the  series 
coefficients  frequently.  This  seems  to  be  a serious  drawback  to 
the  time-series  approach  due  to  the  canputational  requirements  of 
reestimation. 

The  Box-Jenkins  method  for  modeling  a time-series  has  three 
steps:  1)  identification,  2)  parameter  estimation,  and  3) 
diagnostic  forecast  checking.  In  the  identification  process  the 
sample  autocorrelation  is  the  prime  indicator  of  the  order  of 
differencing  required  to  achieve  a stationary  series.  It 
measures  the  correlation  between  K periods  of  the  time-series  art- 
is  computed  by 

n-k 

Dh  - <Vk  - 

-t"'-, 

I (2t  - 

Y » Z - Z = average  deviation  ■' 

T = mean  of  series 
n = number  of  observations  on  series 
k * nunber  of  periods  between  observations 
Z * observation  value  (MW)  from  series 
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Autocorrelations  for  the  new  series  were  examined  to  determine 
which  ARIMA  model  should  tie  initially  tested.  Characteristic 
graphs  of  autocorrelations  (Mabert,  1975)  were  conpared  with 
experimental  results  to  estimate  the  series  order. 

Partial  autocorrelation  is  then  used  to  further  refine  the 
identification  of  a model.  It  measures  the  strength  of  the 
relationship  between  time  periods  in^a  ^series. 


" *kk  >1,2. ..k-1 

During  the  next  step  of  the  Box-Jenkins  method,  estimates  of  the 
coefficients  are  determined  by  the  method  of  least  squares.  For 
a second  order  autoregressive  model 

= autoregressive  coefficient,  i=l  to  order  of  model 

7 = observed  series  value  at  time  t 

t 

Z = forecast  series  value  at  time  t 
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is  minimized.  The  canputational  procedure  involved  uses  the 
identified  autoregressive  model  to  forecast  the  load  Z for  each 
known  observation  in  the  time-series  Z(t).  The  non-linear 
routine  then  calculates  the  forecast  error  at  each  point  and 
selects  the  set  of  which  minimize  this  error.  Finally  the 
model  is  checked  for  adequacy  of  fit  by  using  a chi-square  test, 
Q. 


Q =■  n 


n = total  number  of  observations 
k = number  of  residual  autocorrelation  values 
r (^)  = residual  autocorrelation  of  a at  lag  k 
If  the  canputed  value  of  Q is  greater  than  tabulated  it 
indicates  that  the  model  should  be  rejected  at  that  confidence 
level.  The  IMSL  computer  program  routine  used  provided  the 
significance  level  of  the  chosen  model  directly. 
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III.  1.3  Time-Series  Forecasts 


The  ccmputer  program  (Whitmore,  1977)  was  used  extensively  to 
fit  a time-series  model  to  CVPS's  load  data.  Working  initially 
with  one  months  hourly  data,  weekends  removed,  the  24th 
difference  was  taken  to  reduce  the  data  to  a stationary  series. 


This  produced  an  autocorrelation  of  uniformly  small 
autocorrelation  values  for  increasing  time  lags.  High 
autocorrelation  was  observed  at  k = 1 and  24.  K = 168  hours,  one 
week,  was  another  expected  high  autocorrelation.  The 
autocorrelation  exhibited  an  exponentially  damped  cosine  pattern 
characteristic  of  an  autoregressive  model  with  no  spiking  to 
suggest  a moving  average  conponent.  Applying  the  principle  of 
parsimony  no  further  modification  was  made  to  the  series.  The 
Box-Jenkins  methodology  led  to  the  identification  of  a second 
order  autoregressive  model  from  the  partial  autocorrelation 
vector  as  shown  in  Fig.  13.  The  differenced  model  form  which 
these  values  fit  is: 

p q 

Vd2(t)  “ 

i-1  i=1 


I 


Fig.  13  Partial  Autocorrelation 
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Reducing  t±iis  to  a second  order  undifferenced  autoregressive 
formulation  yields: 


2 (t) 


2 


h (t)  + y Zt.,  (t) 


Canputationally  the  algorithm  iterates  first  on  the  moving 
average  parameters,  then  on  the  autoregressive  parameters  and 
finally  on  the  differences.  Moving  average  and  difference 
operation  were  fixed  to  be  zero.  Hence  a 2nd  order  linear  system 
gives  estimates  of  from  Z=  £ff  Z . Then  a Gaussian  White  Noise 
Test  is  performed  on  the  forecasting  errors  z".  If  these  errors 
pass  the  specified  model  significance  the  model  is  accepted, 
otherwise,  the  program  loops  and  reestimates  0\.  Results  of  a 
typical  forecasting  run  at  the  .02  confidence  level  are  listed  in 
Table  4. 


Of  the  numerious  trial  runs  performed  none  yielded 


TIME-SERIES  FORECAST 


uu 


Constant  = 1 . 3SS 

Autoregressive  Parameters  = .2165  i .4073 


HOUR 

FORECASTS 

FORECASTED 

SERIES 

MAX.  FORECAST 
ERROR  AT  .02 
CONFIDENCE  LEVEL 

1 

12.23 

-27.03 

20.35 

2 

9.83 

-4.67 

20.82 

5 

8.52 

-10.37 

22.78 

4 

7.24 

10.34 

23.10 

5 

5.43 

-14.36 

23.55 

6 

5.73 

.91 

23.69 

7 

5.25 

-16.89 

23.31 

8 

4.36 

-5.32 

23.85 

9 

4.58 

-10.07 

23.90 

10 

4.36 

7.64 

23.91 

1 1 

4. 19 

5.33 

23.93 

12 

4.07 

-10.90 

23.93 

13 

3.98 

5.35 

23.94 

14 

3.91 

-15.89 

23.94 

15 

3.35 

-11.61 

23.94 

16 

3.81 

-22.33 

23.94 

17 

3.78 

-12.01 

23.94 

18 

3.61 

-30.70 

23.94 

19 

3.74 

-3.27 

23.94 

20 

3.73 

-22.30 

23.94 

21 

3.72 

4.42 

23.94 

22 

3.71 

-9.73 

23.94 

23 

3.71 

-7.65 

23.94 

24 

3.70 

5.98 

23.94 

TABLE  4 TIME-SERIES  FORECAST 


L. 
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satisfactory  predictions  for  more  than  four  to  six  hours  in  the 
future  as  can  be  inferred  fran  Table  4.  The  time-series  works 
best  for  short  lead  periods,  and  this  inherent  limitation  has 
been  observed  by  others  (Anderson,  1976) . In  this  case  it  is 
probably  due  to  weather  affects  which  are  not  included  in  the 
time-series  formulation. 


III.  1.4  Frequency  Domain 


It  is  interesting  to  observe  that  due  to  the  periodic  nature 
of  the  basic  system  load  it  could  be  modeled  by  a harmonic 
series.  In  the  frequency  domain  the  load  exhibits  peaks  at 
multiples  of  24  hours  and  168  hours  (one  week) . A Fourier 
analysis  would  give: 


^ (t) 


sin(iwt)  + cosdwt) 


with  w ■ 27^/168  for  the  fundamented  (weekly)  frequency.  The 
problem  then  is  to  estimate  the  coefficients  for  a set  of 
harmonics,  less  than  84,  the  Nyquist  Limit,  so  that  the  periodic 
function  is  adequately  approximated.  Filtering  techniques  have 
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utilized  many  approaches,  but  exponential  smoothing  and  Kalman 
filtering  are  the  most  common  (Gupta,  1972;  Sharma,  1974) . Other 
than  providing  a more  elegant  mathenatical  approach  the 
variations  on  a harmonic  series  model  do  not  claim  to  produce 
average  accuracies  below  4%  or  5%.  It  is  also  apparent  that  the 
operating  branch  of  the  power  conpanies  are  reluctant  to  use  a 
more  conplicated  forecasting  method  than  the  basic  operator's 
prediction.  Other  than  verifying  the  frequency  spectrum  of 
Central  Vermont  Public  Service  Corporation's  load  through 
autocorrelation.  Section  III.  1.2,  no  further  investigation  into 
Fourier  analysis  and  modeling  was  conducted. 

III.  1.5  Linear  Models 

The  operator ' s forecasting  method  determined  a basic  load 
shape  from  years  of  experience.  Often  merely  the  last  day's  or  a 
similar  prior  day's  load  curve  constituted  the  basic  load  shape 
model.  If  the  weather  remained  constant  then  the  basic  load 
shape  became  the  predicted  load.  If  forecasted  weather  differed 
from  a prior  similar  day  then  the  basic  load  shape  had  to  be 
modified.  Operator  forecasting  is  an  art  without  written 
methodology. 


f 
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The  operator's  basic  procedure  generated  a linear  nxxiel  of 
the  form 

Z (t)  - Z (t)  + B (w) 

Z(t)  = Predicted  load 

Z(t)  = hourly  component  from  prior  similar  day 

B(w)  = correction  factor  function  of  weather 
It  is  then  reasonable  to  postulate  a model  vAiich  shifts  Z (t) 
by  an  amount  inversely  proportional  to  the  difference  in 
temperature.  But  what  is  meant  by  the  tanperature  difference? 
Average  temperature  difference  might  be  a better  measure,  however 
small.  The  weather  service  can  not  provide  as  accurate  a short 
term  forecast  of  average  temperature  as  it  can  of  maximum  and 
minimum  temperature. 

The  constraint  of  available  forecasted  tonperature 
information  indeed  iirpacts  on  the  model  selected.  If  a 
tanperature  model  is  chosen  to  follow  the  normal  pattern  of 
morning  low  followed  by  an  afternoon  high  and  the  forecast  is 
also  for  a morning  low  afternoon  high,  then  the  total  model  could 
logically  be  decomposed  into  two  parts.  The  morning  segment  from 


48 


(0001  - 1200  hours)  dependent  on  morning  low  temperatures,  and 
the  afternoon  segment  dependent  on  the  high  forecasted 
temperature. 

Initially  then  a linear  model  might  be: 


(t)  « Z’  (t)  + 
Zg  (t)  - Z«  (t)  + 


4 


z'  = prior  similar  days  load  for  hour  t 

a, = (base  day  AM  low  tanp.  - forecasted  low  tanp.) 

Bj=  (base  day  PM  high  temp.  - forecasted  high  temp.) 

This  formultion  implicitly  assumes  that  the  load  will  change  at  a 
constant  rate  or  gain  factor  of  one  megawatt  per  degree 
Fahrenheit  for  edl  hours  of  the  day  and  year  Fig.  12,  Section  II. 

2.4. 

Using  the  two  years  of  available  data,  results  of  trial 
simulations  produced  average  errors  ranging  frcm  5 - 10%.  Actual 
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maximum  and  minimum  temperatures  were  used  for  these  simulations, 
and  the  results  were  encouraging  when  considered  in  the  context 
of  simplicity  of  procedure.  On  the  basis  of  total  daily  load, 
these  forecasts  were  equivalent  to  those  obtained  by  more 
sophisticated  techniques  in  commerical  practice  (Bechard,  1978). 

To  refine  the  linear  model  into  a dynamic  one  and  at  the  same 
time  make  the  gain  constants  more  realistic,  gain  factors  were 
computed  for  each  hour  from  the  previous  two  similar  days.  This 
formulation  computed: 


( 


G 


i 


A MW  (t^) 

K + [a  Temp  (t^^) 


/ 

G;=  sensitivity  of  load  with  tanperature 
= change  in  system  load  at  hour  t 
T = change  in  temperature  at  hour  t 
K = unity  constant 

Provision  was  made  to  prevent  the  denominator  from  becoming  zero 
by  adding  in  a small  constant.  The  linear  model  equation  then 
became: 
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Z (t)  - Z'  (t)  + (t)  B (t) 

Again  tests  were  tun  on  Central  Vermont  data.  In  many  cases  the 
hourly  results  were  better  than  had  been  observed  with  Gj^=  1. 
However  in  a significant  number  of  samples  a larger  forecast 
error  was  produced.  These  large  errors  corresponded  to  large 
(i.e.  >3)  values  of  G.^.  The  absolute  value  of  was  then 
constrained  less  than  or  equal  to  3 with  a marked  improvement  in 
average  results.  Further  improvements  were  noted  experimentally 
as  Gi  was  reduced  to  a limit  of  1.5.  Beyond  this  value 
prediction  accuracy  again  fell  off.  Experimentally  the  limiting 
value  of  1.5  was  found  to  optimize  forecasting  accuracy  with  this 
cdgorithm  for  the  two  years  of  data  examined. 

However  merely  constraining  the  gain  did  not  insure  that  the  Gj^ 
which  had  not  been  limited  were  optimized  in  a systematic 
fashion.  Recall  that  G^  was  merely  the  ratio  of  change  in  system 
loe»d  to  change  in  tanperature  at  time  t on  two  previous  similar 
days.  While  this  forecasting  procedure  produced  usable  results 
as  the  sample  figures  and  tabulations  of  Appendix  3 illustrate, 
for  non-winter  months,  it  lacked  a formalized  statistical  basis. 
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III.  1.6  Regression  Prediction 

A standard  tool  used  to  make  forecasts  from  statistical  data 
is  regression  analysis.  Straightforward  least  squares  fit  to 
historical  data  merely  reproduce  in  some  manner  the  load  function 
or  curve  shape  vrtiich  is  already  known.  Least  squares 
approximations  alone,  as  a load  forecasting  method,  produce  8%  - 
10%  average  errors  (Galiani,  1975).  A closer  examination  of  the 
linear  regression  methodology  from  SPSS  as  discussed  earlier, 
section  II  2.2,  revealed  two  possible  alternative  regression 
formulations.  The  first  possibility  is  to  perform  the  regression 
analysis  on  each  hour  of  the  prior  thirty  days  producing 
twenty-four  regression  equations; 


(tj^)  - A (t^)  + B (tj^)  Temp  (t^)  i-1,2...24 


Then  the  daily  load  is  merely  the  surmation  of  ).  In 
addition  to  the  necessity  of  reinitializing  constants  A(^  ) and 
B(l^)  frequently,  this  approach  suffered  from  the  practical 
limit,  mentioned  earlier,  that  hourly  tentperature  is  unknown. 
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This  equation  was  modified  to  fit  available  temperature  forecasts 
as  was  done  before.  The  forecasting  equations  then  became: 


'Z  (tj^)  - A (tj^)  + B (tj_)  (T)  1-1, 2. ..12 

(tj)  - A (tj)  + B (tj)  fj  (T)  j-13,14...24 

Z(t  ) = forecasted  load  for  hour  i] 

f^(T)  = prior  day  min.  tanp  - forecasted  min.  temp. 

“ prior  day  max.  temp.  - forecasted  max.  temp. 

Recalling  that  the  basic  linear  regression  equation  has  the  form 

y's  A + BX 

Y'  is  the  estimated  value  of  the  dependent  variable  Y.  The 
constant  A represents  the  Y axis  intercept  of  a best  fit  line 
with  slope  B.  It  would  appear  that  B is  a measure  of  the  gain 
analogous  to  of  the  prior  linear  formulation  section  III. 

1.5.  The  variable  B for  each  hour  measures  the  change  in  Y or 
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system  load  with  a change  of  one  unit  in  X,  the  temperature. 
Implicitly  the  sum  of  squared  residuals,  {^1y-Y'|^  ),  has  been 
minimized  within  the  SPSS  routine  vrtiich  generated  A and  B 
from: 


» Y - B X 


I 


a - I)  (I  - Y) 


- x)2 


where  X and  Y are  mean  values.  How  good  would  a prediction  then 
be  if  maximum  and  minimum  temperatures  were  used  instead  of 
hourly  values?  Table  5 shows  daily  regression  results  for 
October  1976  and  October  1977.  This  entirely  SPSS  method  proved 
to  be  a consistently  poorer  load  predictor  than  did  the  method  of 
Section  III.  1.5 


III.  1.7  Regression-Similar  Day  Method 

A third  alternative  to  the  basic  regression  approach  is  to 
combine  the  regression  slope  constants  B with  a prediction 
routine  based  on  the  previous  similar  day.  In  essence,  the  B 
provide  a best  fit  slope  for  the  ith  hour  prediction  from  the 
prior  thirty  day's  observations. 


COMPARISON  OF  FORECAST 
METHODS 


OCTOBER  1976  AND  OCTOBER  1977 


Average  Errors  1976  Average  Errors  1977 


Date 

Constrained 

SPSS  Gains 

Entirely 

Constrained 

Entirely 

Gains 

Hybrid 

SPSS 

Gains 

SPSS 

1 

2.77 

2.58 

• 

3.4 

16.25 

2 

5.62 

5.29 

20.07 

5.57 

34.17 

3 

8.69 

7.81 

37. 18 

4.91 

4.61 

4 

4.64 

4.67 

10.41 

3.68 

3.52 

5 

4.76 

5. 12 

5.43 

3.  19 

5.12 

6 

2.86 

3.87 

3.39 

3.73 

5.07 

7 

2.2 

2.4 

4. 16 

3.65 

5.75 

3 

3.71 

3.92 

2.11 

8.68 

9.56 

9 

5.22 

5.12 

15.94 

13.32 

16.25 

10 

6.75 

6.39 

26.33 

4.1 

4.3 

1 1 

4.  18 

5.49 

6.34 

6.22 

3.35 

12 

4.71 

3.34 

4.01 

2.02 

— 

13 

3.51 

3.39 

4.4 

3.53 

3.78 

14 

4.05 

2.97 

3.61 

5.54 

4.51 

15 

3.74 

3.59 

2.85 

8.11 

14.11 

16 

7.5 

8.2 

— 

9.95 

25.45 

17 

9.03 

9.83 

37.94 

4.89 

7.26 

18 

3.68 

2.85 

6.27 

5.76 

4.  15 

19 

3.67 

2.82 

3.3 

3.51 

2.93 

20 

5.17 

5.39 

4.29 

3. 16 

2.77 

21 

3.17 

2.69 

4.67 

6.19 

5.2 

22 

4.05 

4.06 

2.23 

5.67 

19.06 

23 

9.98 

9.40 

12.33 

7.41 

33.84 

24 

8.94 

8.39 

— 

4.79 

3.44 

25 

3.96 

4.66 

6.04 

4.77 

3.97 

26 

4.21 

3.09 

3.72 

3.89 

4.21 

27 

4.2 

4.16 

4.6 

3.54 

4.15 

28 

3.25 

3.30 

2.58 

3.95 

5.75 

29 

4.4 

3.2 

4.49 

30 

4.83 

4.93 

17.35 

31 

8.01 

7.26 

30.6 

MEAN 

5.01 

4.84 

10.24 

5.26 

9.35 

STD 

2.06 

2.09 

10.7 

2.4 

9. 14 

DEV 


TABLE  5 Forecasts 
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Since  the  B 's  are  the  least  squares  gains;  their  use  negates 
the  need  to  constrain  the  gain  term  to  be  less  than  1.5.  In 
actuality  their  range  of  values  is  itself  limited  as  can  be  seen 
in  Table  4.  This  model  is  then; 


z (tj)  - 2 (t^)  + 

(T) 

i=1,2...12 

z (tj)  = Z (tj)  + 

(T) 

j=13,14...24 

G » B.  hourly  regression  slope 

38  i 

Canparisons  between  this  hybrid  method  and  that  of  Section 
III.  1.5  were  performed  for  several  days  during  the  peak  months 
in  1976  and  1977  (Appendix  4) . Results  were  consistently 
improved,  approximately  1%,  using  the  hybrid  forecast  algorithm 
for  randomly  selected  days  throughout  both  years.  This  forecast 
algorithm  combines  the  statistical  coefficient  B produced  from 
regression  analysis  and  a functional  weighting  of  the  tanperature 
variable  with  the  prior  similar  day's  load  curve. 

In  reality  the  predicted  load  is  a function  of  a function. 

It  can  be  viewed  as  a function  of  two  variables.  One  variable  is 
a function  of  the  general  load  shape  for  that  day  of  the  week  and 
the  other  variable  is  a function  of  temperature.  Hence,  a model 
can  be  no  better  than  the  basic  function  chosen  to  represent  the 
load  curve  shape.  Regression,  time-series,  sum  of  sinusoids  all 
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approach  the  basic  function  through  a least  squares  averaging 
approximation  which  produces  comparable  results  with  considerably 
greater  conputational  requirements.  It  can  be  legitimately 
argued  that  the  forecast  improvonent  gained  with  the  hybrid  over 
that  of  Section  III.  1.5,  is  not  great  enough  to  offset  the  added 
program  complexity.  Utilities  have  demonstrated  a strong  desire 
to  use  the  simplest  routine  vAiich  will  give  acceptable  results. 

Finally  judging  the  three  methods  on  the  power  industry's 
standard  of  total  required  energy  versus  total  predicted  energy 
each  day,  all  three  methods  are  again  similar  and  equivalent  to 
current  industry  results  (Bechard,  1978). 

III.  2.0  Computer  Algorithm 

The  forecasting  program  written  during  this  study  was 
designed  in  two  independent  parts.  Part  one  is  a general 
interactive  mainline  routine  which  interrogates  the  system 
operator's  console.  This  routine  then  generates  Julian  file 
names  of  similar  days  to  pass  to  a prediction  subroutine.  The 
prediction  subroutines  themselves  constitute  the  second  part  of 
the  forecasting  program.  This  arrangement  allowed  changes  in  the 
prediction  parameters  to  be  readily  made  and  at  the  same  time  be 
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transparent  to  a syston  operator.  Appendix  2 is  a source  code 
listing  of  the  forecasting  program  which  has  been  developed. 

III.  2.1  Similar  Day  Selection 

Conceptually  the  selection  of  a similar  day  is  trivial, 
although  the  similar  day  load  shape  plays  such  a dominant  part  in 
the  forecast  that  several  different  combinations  were  evaluated. 

A best  fit  was  found  to  utilize: 


Day  to  be 

Similar  Day  for  Curve 

Similar  Day  for 

Predicted 

Shape 

Hourly  Gain 

Sunday  (A) 

Prior  Sunday  (A-7) 

Prior  Saturday 

(A-8) 

Monday (B) 

Prior  Monday  (B-7) 

Prior  Friday 

(B-3) 

Tuesday (C) 

Prior  Friday  (C-4) 

Prior  Thursday 

(C-5) 

Wednesday  to 

Prior  Day  (D-1) 

Day  before  Prior 

Friday  (D) 

Day(D-2) 

Saturday (E) 

Prior  Saturday  (E-7) 

Prior  Sunday 

(E-6) 

9 

( 


TABLE  6 Similar  Days 
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Logically  the  prior  similar  day  is  the  most  recent  day  with  the 
same  characteristic  shape  as  the  day  to  be  forecast. 

III.  2.2  Forecast 

Canputationally  the  forecast  is  generated  by  treating  each 
day  as  a 24  component  vector.  This  is  based  on  the  prior  similar 
day  and  modified  by  a scale  factor  or  Gain  operating  on  the 
difference  between  high/low  temperature  on  the  prior  day  and  the 
forecast  temperatures.  The  results  remain  in  a vector  format 
with  each  element  corresponding  to  an  hour.  Sample  tabular 
output  is  in  Appendix  3. 
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CHAPTER  IV 


LOAD  LEVELING 


IV.  1.0  Introduction 

Historical  load  factors  of  Central  Vermont  Public  Service 
Corp.  are  examined.  Impacts  of  deferred  electric  heating  on 
system  load  factors  are  considered  both  with  and  without  load 
forecasts. 


( 


IV.  1.1  Load  Factors  Central  Vermont  Public  Service  Corp. 

The  term  load  factor  only  has  meaning  when  the  time  period 
involved  is  specified.  For  example,  the  daily  load  factor  on  the 
peak  day  in  1976,  (23  January)  is  computed  as; 

Z \ 8132.38 

“ (386.9)  X 24  ' 87. 58^ 

Daily  load  factors  vary  considerably  throughout  the  year 
(Appendix  5) . The  yearly  load  factor  computation  replaces  24 
hours  by  8,760  hours  in  both  the  nunerator  and  denominator,  but 


max 
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the  peak  hourly  term  in  the  denominator  does  not  change. 

IV.  2.0  Electric  Space  Heating 

Mean  hourly  residential  electric  space  heating  consumed 
approximately  3.0  KW/Hr  during  January  1975  (Vermont  Load  Study, 
Chapter  5,  1977) , hence  daily  average  energy  suj^Jlied  to  each 
customer  is  72  KWH.  Therefore,  this  amount  of  average  energy 
must  be  supplied  to  a TES  customer  during  eight  charging  hours 
selected  as  off  peak  times  from  the  forecasted  load.  Assuming 
each  unit  is  charged  at  a rate  of  9 KW/hour  for  the  eight  hours 
it  is  turned  on,  then  one  thousand  customers  will  produce  a 9 MW 
change  in  the  system  load  curves  during  periods  of  TES  operation. 
The  ratio  of  3:1  used  here  is  more  conservative  than  the  industry 
charging  rate  of  2.5. 

It  is  assumed  that  ripple  controlled  TES  units  will 
predominately  be  installed  in  new  dwellings  as  opposed  to 
retrofitting  existing  homes  due  to  their  lower  installation  cost 


in  new  construction. 
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IV.  3.0  Load  Shaping 

For  a given  number,  N in  thousands,  of  new  controlled  TES 
customers  the  load  forecast  curve  of  Section  III.  1.6  or  III.  1.7 
can  first  be  reduced  N x 3 W at  each  of  the  twenty- four  hours. 
This  is  necessary  to  account  for  the  impact  of  these  customers' 
non-deferrable  electric  load.  Then  each  hour  v^ich  is  forecast 
some  arbitary  percentage  below  the  predicted  peak  is  available 
for  operation  of  TES  units.  This  percentage,  assumed  for  this 
study  as  95%,  will  be  determined  by  management  based  upon 
improvement  desired  in  load  factor  and  confidence  in  the 
forecast. 

For  the  purpose  of  this  study  an  hour  is  available  for  TES 
charging  if  the  forecasted  load  is  less  than  or  equal  to  95%  of 
the  reduced  peak.  Then  at  each  available  hour,  Nl,  units  will  be 
operated  to  increase  system  load  to  the  95%  of  the  predicted  peak 
level.  For  the  peak  day  in  1977  a simulation  was  run  with  2,000 
TES  customers.  The  predicted  peak  is  reduced  by  6 MW  to  386.4 
MW,  (actual  peak  was  384. 2rw) , and  twenty  hours  of  predicted  load 
fell  below  the  95%  level.  Now  instead  of  only  eight  to  ten  hours 
availeible  for  nighttime  TES  charging  the  available  charging 
window  has  been  doubled.  With  4,000  customers  18  hours  are 
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available  Fig.  14,  and  daily  load  factor  has  been  improved  to 
93%. 

The  yearly  load  factor  was  only  slightly  improved  to  61%  due 
to  the  summer -winter  load  differential.  Additional  improvements 
to  the  yearly  load  factor  are  possible  with  larger  N,  although 
TES  systems  capable  of  storing  energy  for  six  months  (J.  H.  U. 
Applied  Physics  Laboratory,  1977)  are  more  appropriate  for  annual 
load  leveling.  On  a daily  basis  Central  Vermont  Public  Service 
Corporation  can  consistently  achieve  load  factors  of  over  90%  by 
combining  forecasting,  TES  scheduling,  and  a deferrable  load  of 
70MWH.  This  corresponds  to  approximately  7,500  average  all 
electric  customers. 
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CHAPTER  V 


FURTHER  STUDY 


V.  1.0  Introduction 

Potentially  promising  areas  for  further  forecasting 
investigation  are  proposed  along  with  possible  improvements  to 
the  hybrid  forecasting  cilgorithm  developed  here. 

V.  1.1  Reaxtmendations  for  Further  Study 

No  attempt  has  been  made  in  this  study  to  determine  the 
optimun  number  of  TES  customers  for  a utility.  There  would  seem 
to  be  a limit  beyond  which  additional  TES  units  would  not  be 
econcTrtically  beneficial  from  the  standpoint  of  daily  load  factor 
improvement.  However,  by  controlling  a significant  percentage  of 
winter  loads  the  increased  system  security  is  certainly  valuable 
to  the  utility,  regardless  of  possible  operating  efficiency 
improvements.  In  fact  a completely  TES  systan  is  viable. 
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V.  1.2  Suimier  Peaking  Affects 

The  assumed  daily  temperature  profile  for  this  study  most 
closely  matches  summer  days  and,  hence,  accuracy  is  higher  in 
that  season.  Conversely  it  would  be  valuable  to  study  the 
accuracy  of  this  forecasting  method  when  applied  to  a summer 
peaking  system. 

V.  1.3  Winter  Peaking  Affects 

Could  the  forecasting  model  be  optimized  for  the  heating 
season?  This  is  the  area  of  greatest  potential  improvement  from 
the  standpoint  of  load  leveling  and  probably  of  the  greatest 
benefit  to  most  utilities. 

V.  1.4  Alternate  Applications  of  Forecast 

Assuming  a basic  forecast  is  available  then  to  what  other 
uses  could  a smaill  utility  apply  it?  With  the  advent  of  smaller 
municipally  owned  utilities  vAiich  control  their  own  limited 
generating  facilities,  what  accuracy  forecast  versus  degree  of 
ccniplexity  will  be  required  for  their  efficient  and  secure 
operation? 
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Buying  and  selling  blocks  of  power  generating  capacity 
between  utilities  is  often  done  on  contracts  lasting  only  a few 
days  to  a month.  As  fossil  energy  costs  continue  to  rise  the 
selling  of  alternate  capacity  may  increase.  Can  the  algorithm- be 
extended  to  provide  a basis  for  these  sales  by  forecasting 
surpluses? 

V.  2.0  Recommended  Model  Improvements 

Other  weather  variables  such  as  wind  speed  and  cloud  cover 
which  were  not  included  in  the  current  forecasting  model  are 
readily  available  in  the  existing  data  files  generated  in  this 
study.  Further  study  incorporating  cloud  and  wind  into  the 
prediction  algorithm  would  probably  improve  accuracy. 

A multivariate  time-series  program  is  now  available  on  the 
University  computer  system.  Improved  accuracies  by  including 
weather  variables  in  the  time-series  formulation  are  probable  and 
should  be  investigated. 
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CHAPTER  VI 


CONCLUSION 


Two  conclusions  can  be  drawn  frcm  this  study.  First,  that 
systan  loads  such  as  those  of  Central  Vermont  Public  Service 
Corporation  can  be  simply  forecast,  using  a first  order  linear 
model.  Average  yearly  forecast  errors  are  approximately  5%.  The 
ccmputation,  data  storage,  and  software  maintenance  requirements 
are  minimal,  so  that  a small  utility  can  utilize  the  algorithm  at 
the  operating  level  without  additional  dedicated  resources. 

Secondly,  operating  economies  \4iich  can  be  derived  from  a 
foreknowledge  of  system  loads  extend  beyond  mere  scheduling  of 
generating  capacity,  but  includes  expanded  potential  for  systan 
security  through  controlled  load  shedding  using  ripple  control. 
Combining  a forecast  with  the  availability  to  defer  loads  (TES) 
and  centralized  control,  inmediately  improves  system  load 
factors.  A utility  with  these  three  capabilities  will  possess 
distinct  operating  advantages  in  the  energy  scarce  future. 
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APPENDIX  1 


ORGANIZATION  OF  DATA 


The 
files  is 


standard  annual  calendar  used  to  name  data 
presented  followed  by  a sample  day' s data. 


Cloud  cover  is  coded: 

0. . . .clear 

2.. .. partly  cloudy 

4.. .. cloudy 

6. .  . .rain 


8. . . .snow 
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The  Number  of  Each  Doy  of  the  Yeor 


259  289 

260  290 

2C1  291 

262  292 

263  293 

264  294 

265  295 

266  296 

267  297 

268  298 

269  299 

270  300 

271  301 

272  302 
TTZ  303 

304 


305  335 

306  336 

3C7  337 

308  338 

309  339 

310  340 

311  341 

312  343 

313  343 


315  345 

316  346 

317  347 

318  348 

319  349 

320  350 

321  351 

322  352 

323  353 

324  354 

325  355 

326  356 

327  357 

328  358 

329  359 

330  360 

331  361 

332  362 

333  363 

334  364 
365 


For  leap  years  the  number  of  the  day  after 
February  28  is  one  greater  than  that  given  in  the  table. 
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APPEMDIX  2 

COMPUTER  PRC3RAM3 

The  programs  listed  are  written  in  Fortran  IV  for  the 
university  of  Vermont  Sigma  5 operating  system.  The  library 
routines  utilized  were  current  on  the  system  as  of  the  date  of 
this  report. 

Pages  74  to  75  contain  a listing  of  a general  mainline 
program  which  interfaces  the  system  operator  with  a particular 
forecasting  program.  This  routine  also  selects  two  prior  similar 
days  on  which  a forecast  is  based. 

Pages  77  to  20  are  a subroutine  which  forecasts  the  hourly 
system  load  twenty-four  hours  in  advance  using  gains  fi.xed  less 
than  1.5. 

Pages  31  to  34  are  another  subroutine  which  forecasts  the 
system  load  for  February  1975,  but  which  employs  manually  entered 
hourly  gain  coefficients  produced  off-line  v;ith  a statistical 
routine . 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

C PROGRAM  TO  INTERFACE  DISPATCHERS  AND  FORECASTING  ROUTINES 

C 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C PROGRAM  MUST  BE  USED  WITH  A PREDICTION  SUBROUTINE 
C 

C FORT  OPTIONS  MUST  BE  BOTH  BC  AND  S 
C 

C LINK  WITH  : LIB. EDS 
C 

DIMENSION  MWARRAY(24,1) ,TEMPARRAY(24, 1) ,JULIANDATE(3) ,D(10) ,A(12) 
DIMENSION  FILE(3) ,DIGET(3) ,DATE(2) ,DATE3(2) ,DATE2(2) ,DATE1(2) 
INTEGER  CHANGE  1 .HOLIDAY, CHANGE2,YEARCHANGE 
IMPLICIT  INTEGER  (A-H) 

COMMON  YEA RCHANGE, CHANGE 1,CHANGE2 
CALL  LIBINIT 
C 

C READ  DATE  TO  BE  FORECAST  FROM  CRT 

WRITE(108,10) 

10  FORMATCWHAT  IS  THE  JULIAN  DATE  TO  BE  FORECAST  ??  I. E.  100-76') 

C 

READ  (108,20)  DATE3(I),I=1,2 
20  FORMAT (2A4) 

WRITE  (108,30) 

C READ  DAY  OF  WEEK  TO  BE  FORECAST  FROM  CRT 

C 

30  FORMATCWHAT  DAY  OF  WEEK  IS  TO  BE  FORECAST??') 

C 

READ  (108,40)DAYOFWK 
40  FORMAT  (A2) 

C 

C CHECK  TO  SEE  IF  FORECAST  IS  FOR  A SPECIAL  HOLIDAY 
WRITE  (108,50) 

50  FORMAT  ('IS  THE  DAY  A SPECIAL  HOLIDAY  OTHER  THAN  SAT  OR  SUN?') 

READ  (108, 60) HOLIDAY 
60  FORMAT  (A4) 

C 

C SET  ERROR  FLAG 
CHANGE 1=2 
C 

C SET  BRANCH  CHANGES  FOR  APPROPRIATE  SIMILIAR  DAY 
C 

IF  (DAY0FWK.EQ.2HSU)  CHANGE 1=-7;CHANGE2=-1 

IF  (DAY0FWK.EQ.2HSA)  CHANGE 1=-7;CHANGE2=+1 
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IF  (DAY0FWK.EQ.2HTU)  CHANGE1=-J4;CHANGE2=-1 

C 

IF  ( DA YOFWK . EQ . 2HMO ) CHANGE  1 =-7 ; CHANGE2=+4 
C 

IF  ((DA YOFWK  .EQ.  2HWE)  .OR.  (DAYOFWK  .EQ.  2HTH)  .OR. 
#(DAYOFWK  .EQ.  2HFR))  CHANGE  1 =- 1 ;CHANGE2r-1 
C 
C 

YEARCHANGEsO 

IF(HOLIDAY  .EQ.  1HY)YEARCHANGE=-1 ;CHANGE1=0;CHANGE2=-1 
C 

C STOP  IF  ERROR  FLAG  HAS  NOT  BEEN  CHANGED  FROM  2 
C 

IF  (CHANGE1  .EQ.  2)  WRITE  (108,70);GO  TO  101 
70  FORMAT  ('INPUT  ERROR  IN  DATE  DAY  OR  HOLIDAY.  JOB  STOPPED!') 

C 

C COMPUTE  NAME  OF  PRIOR  SIMILIAR  DAY  FILE 
C 

DATEKI  )=DATE3(1 ) 

DATE1(2)=DATE3(2) 

CALL  BASEDATE  ( DATE 3, CHANGE  1 , YEA RCHANGE) 

C 

C COMPUT  NAME  OF  EARLIEST  SIMILIAR  DAY  FILE  FOR  GAINS 
C 

DATE2(1 )=DATE3(1 ) 

DATE2(2)=DATE3(2) 

CALL  BASEDATE  (DATE3,CHANGE2,YEARCHANGE) 

C OUTPUT  TO  CRT  NAMES  OF  3 FILES  USED  (OPTIONAL) 

C 

WRITE  (108,9998)  DATE1 ,DATE2,DATE3 
C 

9998  FORMAT  (X,  'DATEU*  ,X,2A4,4X,  'DATE2='  ,X,2A4,4X,  'DATE3=’  ,2A4) 
CALL  PREDICT  (DATE  1 , DATE2, DATE3) 

101  CALL  EXIT 

STOP  'NORMAL  ENDING’ 

END 

C 

C 

C SUBROUTINE  TO  COMPUTE  EARLIER  DAY  FILE  NAME 

C 

SUBROUTINE  BASEDATE  ( DATE, CHANGE, YEARCHANGE) 

INTEGER  DATE(2) , TEMP, CHANGE, YEARCHANGE 
C CONVERT  NAME  TO  NUMERIC 
C 

t 


s 

LI,1 

1 

s 

LW,3 

•DATE,  1 

s 

LI, 2 

X'FO’ 

s 

SLD,2 

-8 

s 

STW,3 

•DATE, 1 

r 
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DECODE(4,9,DATE(1 ))DATE(1 ) 

DATE ( 1 ) =D ATE ( 1 ) +CHANGE 
C 

C CHECK  TO  SEE  IF  YEAR  COMPUTED  IS  CORRECT 
C 

IF  (DATE(I)  .LE.  0)  DATE( 1)=365+CHANGE ;YEARCHANGE=-1 
IF  (DATE(I)  .GE.  367)  STOP  'JULIAN  DATE  ERROR* 

C 


C 

C 


C 

C 

C 

102 

S 

S 

S 

S 

S 


9 

90 

99 

999 

997 

998 
9999 


TEMPrDATEd  ) 

MAKE  SURE  PROPER  NUMBER  OF  ZERO'S  ARE  IN  NAME 

IF(DATE(1).LE.9)ENCODE(n,997,DATE(1))TEMP;GO  TO  102 
IF(DATE(1 ) .LE.99)ENC0DE(4,998,DATE(1 ))TEMP 

CONVERT  BACK  TO  ALPHABETIC  NAME 

IF  (DATE(I)  .GE.  100 )ENCODE ( 4 , 999 , DATE ( 1 ) )TEMP 
LI,1  2 

LW,3  »DATE,1 

LI, 2 X'FOFO' 

SLD,2  -16 

STW,3  . *DATE,1 
DECODE (4,9, DATE ( 2 ) ) DATE ( 2 ) 

TEMP=DATE(2) 

TEMPrTEMP+YEARCHANGE 
ENCODE (4,90, DATE ( 2 ) ) TEMP 
FORMAT (14) 

FORMAT (12) 

F0RMAT(X,A4) 

F0RMAT(I3, '-') 

FORMAT  COO'  ,11, '-') 

FORMATCO'  ,12,'-') 

F0RMAT(X,2A4) 

RETURN 

END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

C LINEAR  PREDICTION  ROUTINE 

C GAINS  COMPUTED  FROM  PREVIOUS  HOUR  OF  PREVIOUS  TWO  DAYS 

C 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C CALLING  PROGRAM  PASSES  THREE  FILE  NAMES  TO  THIS  SUBROUTINE 
C ALL  FILES  NAMED  MUST  BE  ON  DISK  OR  THE  OPENF  CALLS  CHANGED 
C TO  ACCESS  TAPES. 

C OUTPUT  CONSISTS  OF  TWO  FILES  'ACTUAL'  AND  'PREDICTED'  LOADS 
C IN  ADDITION  TO  LP  4 CRT  TABULATIONS. 

C 

C DCB#3  IS  THE  FID  TO  BE  PREDICTED 

C DCBtf2  IS  THE  FID  FOR  PRIOR  SIMILIAR  DAY 

C DCB#2  IS  THE  FID  FOR  EARLIEST  DAY  USED  IN  GAIN  PREDICTION 

C 

C SUBROUTINE  MUST  BE  USED  WITH  PROGRAM  'MAINLINE'  AND 
C LINKED  WITH  : LIB. EDS 

C FORT  OPTIONS  MUST  BE  BC,S 
C 
C 

SUBROUTINE  PREDICT  ( D ATE  1 , DATE 2, DATES) 

DIMENSION  DAY(24,2) ,FUTURE(24,2) ,PREDICTI0N(24) ,DIFFERENCE(24) , 
* PERCENT ( 24 ), GAIN ( 24 ) , DA YONE (24,2), DENOMATOR ( 24 ) 

REAL  PMFACT 

INTEGER  H1,H2,H3 

CALL  OPENF  ( 1 ,DATE3, 2,0,0, 1 , 1 ,7,0,0,0) 

CALL  OPENF  (2,DATE2,2,0,0, 1 , 1 ,7,0,0,0) 

CALL  OPENF  (3,DATE1 , 2, 0, 0, 1 , 1 , 7, 0, 0, 0) 

C 

C READ  MW  & TEMPERATURE  BY  HOUR  INTO  3 ARRAYS  FOR  3 DAYS  USED 
C 

DO  40  1=1,24 

READ  (1,30)  (DAY0NE(I,J),J=1,2) 

READ  (2,30)(  DAY(I,J),J=1,2) 

READ  (3,30)  (FUTURE(I,J),J=1,2) 

30  F0RMAT(X,T35,F5.1,T50,G) 

40  CONTINUE 

C 

C READ  AND  OUTPUT  DAY  AND  DAY  INFO  ON  DAY  TO  BE  PREDICTED 
C 

REWIND(3) 

READ  (3,31)  PTE31 ,PTE32,PTE33,PY3 
FORMAT  (3A4,T20,A4) 

WRITE  (6,32)  PY3,PTE31 ,PTE32,PTE33 
WRITE  (103,32)  PY3,PTE31 ,PTE32,PTE33 


31 
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32  FORMAT  (1H1, 'PREDICTION  FOR ' ,2X , A4 ,X, 3AU . 2X ' 

C 

C READ  AND  OUTPUT  DATE  i DAY  INFO  ON  PRIOR  SIMILIAR  DAY 
C 

REWIND(2) 

READ  (2,31)  PTE21 ,PTE22,PTE23,PY2 
WRITE  (6,33)  PY2,PTE21 ,PTE22,PTE23 
WRITE  (108,33)  PY2,PTE21 ,PTE22,PTE23 

33  FORMAT(X, 'PREDICTION  USES' ,2X,A4,X, 3AU,X, ' AS  BASE  DAY.',/) 

C 

C READ  AND  OUTPUT  DATE  & DAY  INFO  ON  EARLIEST  DAY  FOR  GAIN  ' 

C 

REWIND(I) 

READ  (1,31)  PTE11 ,PTE12,PTE13,PY1 
WRITE  (6,34)  PY1 ,PTE11 ,PTE12,PTE13 
WRITE  (108,34)PY1 ,PTE11 ,PTE12,PTE13 

34  FORMAT  (X,'GAIN  BASED  ON' ,2X,A4,X, 3A4,/) 

C 

C CHECK  FILES  FOR  HOLIDAY 

C 

REWIND  (1) 

READ  (1,35)  HI 
REWIND  (2) 

READ  (2,35)  H2 
REWIND  (3) 

READ  (3,35)  H3 

35  FORMAT  (X,T12,A4) 

IF  ((H1.EQ.1HH).OR.(H2.EQ.1HH).OR,(H3.EQ.1HH))WRITE  (6,37); 
#WRITE  (108,37) 

37  FORMAT  (X, '»•»»»•  WARNING  ONE  OF  THE  DAYS  USED  IN  THIS  FORECAST 
#WAS  A HOLIDAY  ••»»•*',/) 

C 

C PRINT  INFO  ON  ANY  GAIN  RESTRICTION  FOR  REF. 

C 

WRITE  (6,38) 

38  FORMAT (X, 'GAIN  IS  FIXED  LESS  THAN  1.5.  ' ,/) 

C 

C FIND  MIN.  i MAX.  TEMPERATURES  FROM  PRIOR  SIMILIAR  DAY 
C 

TEMPMAX=-20.0 
TEMPMIN =400.0 
C 

DO  60  1=1,24 
C 

IF(DAY(I,2)  .GT.  TEMPMAX)  TEMPMAX=DAY(I , 2) 

IF(DAY(I,2)  .LT.  TEMPMIN)  TEMPMIN=DAY(I,2) 

C 

60  CONTINUE 

C 

C INPUT  FROM  CRT  FORECASTED  MIN.  & MAX.  TEMPERATURES 


c 
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WRITE  (108,98) 

98  FORMAT  (X.’WHAT  IS  THE  MIN.  FORECASTED  TEMPERATURE?') 

READ  (108,97)  FUTUREMIN 
97  FORMAT (G) 

C 

WRITE  (108,96) 

96  FORMAT(X, 'WHAT  IS  THE  MAX.  FORECASTED  TEMPERATURE?') 

READ  (108,97)  FUTUREMAX 
C 

C COMPUTE  HOURLY  GAINS  G(I) 

C 

DO  85  1=1,24 

DENOMATOR (I)=DAY(I,2)-DAY0NE(I,2) 

IF  (DENOMATOR(I)  .EQ.  0.0)  DENOMATOR ( I )= 10.0 
C 

C RESTRICT  GAINS  LESS  THAN  1.5 
C 

GA IN ( I ) = ABS ( ( DA Y ( 1 , 1 ) -DA YONE (1,1)) /DENOMATOR ( I ) ) 

IF  (GAIN(I)  .GE.  1.5)  GAIN(I)=1.0 

85  CONTINUE 
C 

C MORNING  PREDICTION 

C 

SCALEAMs (TEMPMIN-FUTUREMIN ) 

C 

C PREVENT  DENOMINATOR  FROM  = 0 

C 

IF  (SCALEAM  .EQ.  0.0)  SCALEAMsI.O 
DO  80  1=1,12 

PREDICTION(I)=(DAY(I, 1))+((GAIN(I))»(SCALEAM)) 

80  CONTINUE 

C 

C AFTERNOON  LOAD  PREDICTIONS 

C 

CALL  OPENF( 4, 'ACTUAL  ' ,4, 0,0, 0,0, 7) 

CALL  OPENF  (5, 'PREDICTED  ' ,4, 0,0, 0,0, 7) 

SCALEPMs ( TEMPMAX-FUTUREMAX ) 

IF  (SCALEPM  .EQ.  0.0)  SCALEPM=1.0 
C 

DO  100  1=13,24 

PREDICTION ( I ) = (DA Y ( I , 1 ) ) f ( (GAIN ( I ) ) • ( SCALEPM) ) 

100  CONTINUE 

TOTAL =0.0 
C 

C OUTPUT  RESULT  HEADINGS  TO  LP  & CRT 
C 

WRITE  (108,95) 

WRITE  (6,95) 

95  FORMAT (X, 'PREDICTION' ,T16, 'ACTUAL  MW' ,T30, 'DIFFERENCE' ,153,'*’. 
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•T63, 'GAIN*  ,/) 

TOTALMWPREsO.O 
TOTALMWFUTsO.O 
DO  110  1=1, 2n 

D1FFERENCE( I )= (PREDICTION( I )-FUTURE(I , 1 ) ) 

PERCENT ( I )= ( ABS ( DIFFERENCE ( I ) ) /FUTURE ( 1 , 1 ) ) * 1 00 
WRITE ( 1 08 , 99 )PREDICTION ( I ) .FUTURE (1,1), DIFFERENCE (I ) , PERCENT ( I ) 
•,GAIN(I) 

WRITE(6,99)PREDICTION(I),FUTURE(1, 1)  .DIFFERENCEd)  .PERCENTd) 
»,GAINd) 

99  F0RMAT(X,F8.4,T16,F8.4,T30,F8.4,T50,F6.2,T60,F8.4) 

TOTALsTOTAL+PERCENT ( I ) 

TOTALMWPRErPREDICTION ( I )+TOTALMWPRE 
TOTALMWFUTsTOTALMWFUT+FUTUREd,  1 ) 

WRITE  (4,90)I,FUTUREd,  1) 

90  FORMAT  (X,I,F6.2) 

WRITE  (5,90)  I, PREDICTION (I) 

110  CONTINUE 

CALL  CL0SEF(4,2) 

CALL  CL0SEF(5,2) 

COMPUTE  & OUTPUT  TOTAL  DAILY  ERROR 

DAYPERCENT=(TOTALMWFUT/TOTALMWPRE)»100 
AVGPERCENTr  TOTAL/24.0 
WRITE  (108,999)  AVGPERCENT.DAYPERCENT 
WRITE  (6,999)  AVGPERCENT.DAYPERCENT 
99  FORMAT  (/,’AVG.  OF  HOURLY  PERCENT  ERRORSs' ,2X,F6.2,//, 

# 'PERCENT  ERROR  FOR  TOTAL  DAILY  POWER=' ,2X,F6.2) 

CLOSE  ALL  FILES 

CALL  CLOSEF  (1,2) 

CALL  CLOSEF  (2,2) 

CALL  CLOSEF  (3,2) 

RETURN 
END 
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ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

C LINEAR  PREDICTION  ROUTINE 

C GAINS  COMPUTED  BY  SPSS  SLOPE  FOR  30  DAY  PERIOD 

C »»*«  fEB  76  *••• 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C CALLING  PROGRAM  PASSES  THREE  FILE  NAMES  TO  SUBROUTINE 

C ALL  FILES  NAMED  MUST  BE  AVAILABLE  ON  DISK  OR  OPENF 

C ROUTINES  CHANGED  TO  ACCESS  TAPES. 

C OUTPUT  CONSISTS  OF  LP  & CRT  TABLES.  TWO  FILES  ARE  ALSO 

C CREATED  CONSISTING  OF  'ACTUAL'  AND  'PREDICTED'  LOADS. 

C DCB#3  IS  THE  FID  TO  BE  PREDICTED 

C DCB#2  IS  THE  FID  FOR  THE  PRIOR  SIMILIAR  DAY 

C DCB#1  IS  FID  FOR  EARLIEST  DAY  USED  IN  GAIN  COMPUTATION 

C 

C SUBROUTINE  MUST  BE  USED  WITH  PROGRAM  'MAINLINE'  AND 

C LINKED  WITH  : LIB. EDS 

C FORT  OPTIONS  MUST  BE  DECLAIRED  BOTH  BC  AND  S. 

C 

C 

SUBROUTINE  PREDICT  (DATE1 ,DATE2,DATE3) 

DIMENSION  DAY(24,2),FUTURE(24,2) ,PREDICTION(24),DIFFERENCE(24) , 
• PERCENT(24) ,GAIN(24),DAY0NE(24,2) ,DEN0MAT0R(24) 

REAL  PMFACT 

INTEGER  H1,H2,H3 

CALL  OPENF  (1 ,DATE3, 2,0,0, 1 , 1 ,7, 0,0,0) 

CALL  OPENF  (2,DATE2,2,0,0, 1 , 1 ,7,0,0,0) 

CALL  OPENF  (3,DATE1 ,2,0,0, 1 , 1 ,7,0,0,0) 

C 

C READ  MW  & TEMP  BY  HOUR  INTO  3 ARRAYS  FOR  3 DAYS  USED 
DO  40  1=1,24 
C 

READ  (1,30)  (DAY0NE(I,J),J=1 ,2) 

READ  (2,30)(  DAY(I,J),Js1,2) 

READ  (3.30)  (FUT0RE(I,J),J=1,2) 

30  F0RMAT(X,T35,F5.1,T50,G) 

40  CONTINUE 

C 

C READ  AND  OUTPUT  DATE  AND  DAY  INFO  ON  DAY  TO  BE  PREDICTED 
REWIND(3) 

READ  (3,31)  PTE31 ,PTE32,PTE33.PY3 

31  FORMAT  (3A4,T20,A4) 

WRITE  (6,32)  PY3,PTE31 ,PTE32,PTE33 
WRITE  (108,32)  PY3,PTE31 ,PTE32,PTE33 

32  FORMAT  (1H1, 'PREDICTION  FOR', 2X,A4,X,3A4,2X, '*»••' ,/) 

C 
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C READ  AND  OUTPUT  DATE  i DAY  INFO  ON  PRIOR  SIMILIAR  DAY 
REWIND(2) 

READ  (2,31)  PTE21 ,PTE22,PTE23,PY2 
WRITE  (6,33)  PY2,PTE21 ,PTE22,PTE23 
WRITE  (108,33)  PY2,PTE21 ,PTE22,PTE23 

33  FORMAT(X, ’PREDICTION  USES’ ,2X,A4,X, 3A4,X, 'AS  BASE  DAY.’,/) 
C 

C READ  AND  OUTPUT  DATE  i DAY  INFO  ON  EARLIEST  DAY  FOR  GAIN 
REWIND(I) 

READ  (1,31)  PTE11 ,PTE12,PTE13,PY1 
WRITE  (6,34)  PY1,PTE11 ,PTE12,PTE13 
WRITE  (108,34)PY1 ,PTE11 ,PTE12,PTE13 

34  FORMAT  (X,’GAIN  BASED  ON’ ,2X,A4,X, 3A4,/) 

C 

C CHECK  EACH  DAY  FOR  HOLIDAY 
REWIND  (1) 

READ  (1,35)  HI 
REWIND  (2) 

READ  (2,35)  H2 
REWIND  (3) 

READ  (3,35)  H3 


35  FORMAT  (X,T12,A4) 

IF  ((H1.EQ.1HH).0R.(H2.EQ.1HH).0R.(H3.EQ.1HH))WRITE  (6,37); 
#WRITE  (108,37) 

37  FORMAT  (X, ’*••*»•  WARNING  ONE  OF  THE  DAYS  USED  IN  THIS  FORECAST 
#WAS  A HOLIDAY  »»*•••’,/) 

C 

WRITE  (6,38) 

38  FORMAT (X, ’GAINS  COMPUTED  BY  SPSS  FOR  MONTH  OF  PREDICTION’,/) 
C 

C FIND  MIN.  & MAX.  TEMPERATURE  FROM  PRIOR  SIMILIAR  DAY 


TEMPMAX=-20.0 

TEMPMIN=400.0 

C 

DO  60  Ir1,24 
C 

IF(DAY(I,2)  .GT.  TEMPMAX)  TEMPMAX=DAY(I,2) 
IF(DAY(I,2)  .LT.  TEMPMIN)  TEMPMIN=DAY(I,2) 

C 

60  CONTINUE 

C 

C INPUT  FROM  CRT  FORECASTED  MIN.  & MAX.  TEMPERATURES 


WRITE  (108,98) 


98 

FORMAT  (X,’WHAT  IS  THE  MIN. 
READ  (108,97)  FUTUREMIN 

FORECASTED 

TEMPERATURE?’ ) 

97 

FORMAT (G) 

WRITE  (108,96) 

96 

FORMAT(X, ’WHAT  IS  THE  MAX. 

FORECASTED 

TEMPERATURE?’ ) 

READ  (108,97)  FUTUREMAX 

o o 
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C SPSS  GAINS  LISTED  HERE  FOR  MONTH  OF  PREDICTION  (FEB  76) 

GAIN(I)  =1.5 
GAIN(2)  =.989 
GAIN(3)  =1.19 
GAIN  (14)  =.56 
GAIN(5)  =.624 
GAIN(6)  =.474 
GAIN(7)  =1.18 
GAIN(8)  =.21 
GAIN(9)  =.62 
GAIN(10)=.24 
GAIN(11)=1.6 
GAIN(12)=1 .08 
GAIN(13)=1. 16 
GAIN(14)=.752 
GAIN(15)=.389 
GAIN(16)=1 .07 
GAIN(17)=.81 
GAIN(18)=1.02 
GAIN(19)=.84 
GAIN(20)=.52 
GAIN(21)=1.03 
GAIN(22)=.196 
GAIN(23)=.995 
GAIN(24)=1. 12 

C MORNING  LOAD  PREDICTIONS 

C 

SCALEAM=(TEMPMIN-FUTUREMIN) 

IF  (SCALEAM  .EQ.  0.0)  SCALEAMsI.O 
DO  30  1=1,12 

PREDICTION ( I )= (DAT ( I , 1 ) ) + ( (GAIN ( I ) ) * ( SCALEAM) ) 

80  CONTINUE 

CALL  OPENF( 4, 'ACTUAL  ’ ,4, 0,0, 0,0, 7) 

CALL  OPENF  (5, 'PREDICTED  • ,4, 0,0, 0,0, 7) 

C AFTERNOON  LOAD  PREDICTIONS 

C 

SCALEPM=(TEMPMAX-FUTUREMAX) 

IF  (SCALEPM  .EQ.  0.0)  SCALEPM=1.0 
C 

DO  100  1=13,24 

PREDICTI0N(I)=(DAY(I,1))+((GAIN(I))»(SCALEPM)) 

100  CONTINUE 

T0TAL=0 . 0 

OUTPUT  RESULT  HEADINGS  TO  LP  & CRT 
WRITE  (108,95) 

WRITE  (6,95) 

95  FORMAT (X, 'PREDICTION' ,Tl6, 'ACTUAL  MW' ,T30 ,' DIFFERENCE' ,T53,'S'  , 
*T63, 'GAIN' ,/) 


no  o o 
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COMPUTE  i ERROR  FOR  EACH  HOUR 

TOTALMWPREsO.O 
T0TALMWFUT=0.0 
DO  110  1=1,24 

DIFFERENCE(I)=(PREDICTION(I)-FUTURE(1, 1 )) 

PERCENT ( I ) = ( ABS (DIFFERENCE ( I ) ) /FUTURE ( 1 , 1 ) ) * 1 00 
OUTPUT  RESULTS  TO  LP.CRT  AND  OUTPUT  FILES 

WRITE(108,99)PREDICTI0N(I) ,FUTURE(I, 1),DIFFERENCE(I),PERCENT(I) 
»,GAIN(I) 

WRITE(6,99)PREDICTION(I)  .FUTUREd,  1)  ,DIFFERENCE(I),PERCENT(I) 
*,GAIN(I) 

99  FORMAT(X,F8.4,T16,F8.4,T30,F8.4,T50,F6.2,T60,F8.4) 

TOTALrTOTAL+PERCENT ( I ) 

T0TALMWPRE=PREDICTI0N(I)+T0TALMWPRE 
TOTALMWFUT=TOTALMWFUT+FUTURE ( 1,1) 

WRITE  (4, 90)1, FUTUREd,  1) 

90  FORMAT  (X,I,F6.2) 

WRITE  (5,90)  I,PREDICTION(I) 

110  CONTINUE 

CALL  CL0SEF(4,2) 

CALL  CLOSEF(5,2) 

C 

C COMPUTE  AND  OUTPUT  TOTAL  DAILY  ERRORS 

DAYPERCENTs ( TOTALMWFUT/TOTALMWPRE ) * 1 00 
AVGPERCENTs  TOTAL/24. 0 
WRITE  (108,999)  AVGPERCENT,DAYPERCENT 
WRITE  (6,999)  AVGPERCENT,DAYPERCENT 
999  FORMAT  (/,'AVG.  OF  HOURLY  PERCENT  ERRORS:' ,2X,F6. 2,//, 

f 'PERCENT  ERROR  FOR  TOTAL  DAILY  POWER:' ,2X,F6.2) 

C 

C CLOSE  ALL  FILES 

CALL  CLOSEF  (1,2) 

CALL  CLOSEF  (2,2) 

CALL  CLOSEF  (3,2) 

RETURN 

END 
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SAMPLE  RESULTS 


This  appendix  contains  graphical  follov 
of  thirty* load  forecasting  simulations  on  1 
dates  used  were  selected  by  an  APL  random  r 
and  include  weekends.  For  all  graphs  the  s 

— 0 — Actual  load  on  forec 

— — Predicted  load 

The  mean  error  for  these  30  days  was  5. 
deviation  of  2.68. 
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PRtOlCTIPN  F‘^R  TUtS  10 

76 

PRtOlCTIPN  DSFS  FWI  6 FrB 

76  AS  BASF  Day. 

Gain 

Based  9n 

ThUK  5 FEd  76 

gain 

IS  FIXED 

less  Than  1.5» 

PRLDICTIPN 

AcTliAl  YIW 

OlFFERENfF 

% 

21». 

1000 

224 .booo 

-6.4000 

2 

.35 

20fe« 

7000 

^OAt^OO^ 

**.8000 

2 

.38 

222* 

3000 

213.3000 

8.9000 

a 

.17 

2Vd» 

‘*000 

217.1000 

-4.7000 

2 

*16 

20a» 

7000 

212.8000 

-A. inoo 

1 

.93 

E28» 

1800 

219.2000 

P.9SO0 

4 

•10 

26y* 

bOOl 

261.8000 

7.8000 

2 

.98 

31/» 

o,aQij 

303.0000 

14.0000 

4 

.•  p2 

33u* 

Rooo 

321 .7000 

B .8000 

2 

.74 

32‘A* 

3000 

30A.8O00 

l9.b0OO 

6 

.40 

320* 

3999 

311.3999 

9.0000 

2 

.SP 

33u« 

7000 

299.8999 

30*8000 

10 

.27 

282  > 

3999 

288.1001 

-5. 7oo2 

1 

.38 

268  . 

7QQQ 

277. 7nrtn 

-9.0000 

3 

. P4 

272. 

booo 

273.3999 

-2.8999 

1 

*0= 

262* 

0000 

26b .3000 

-3.3000 

1 

.?4 

268. 

3000 

2S2.1001 

-13.8000 

4 

. 39 

26y. 

bAo'* 

290.3999 

-20*7595 

7 

*15 

298. 

lOOl 

289.2000 

8.9001 

3 

*08 

23‘^» 

0000 

281 .booo 

12. booo 

4 

.44 

281. 

3999 

274.0000 

7.3999 

2 

•70 

264. 

booo 

250*9000 

13.6000 

5 

. 42 

25/. 

iJDOi 

2‘t3.3Q00__ 

13. 8001 

5 

.67 

21b. 

7000 

20*+. booo 

12*1000 

5 

.91 

AVG. 

8F  Hourly 

Emergent  err'^rr 

= 3.93 

PERC.ENT  EpRBR 

FBR  total  daily 

P9WtR»  98*28 

BBUVAIttBLLCOPY 
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P«tDlCTI'3N  FSP  wru  1 i<  F£3  7© 

PHLOICTI^^N  U3FS  TU£S  17  Ff6  76  AS  BASE  ^AY. 


GAIN  Based  0N  M3N  is  FeB  76 


Gain  is  t^jxED 

LE3S  than 

1 .s. 

PHLDICT 1 bn 

actual 

OIPFEPENCE 

% 

la/ 

1667 

201 

•700C 

-14. 9333 

7.21 

isy 

0000 

189 

•7000 

-.7300 

.37 

I7d 

3B71 

1 90 

•7000 

-17.3429 

?.09 

l7b 

4Sb7 

181 

*4030 

-5.9143 

3.26 

i6y 

7167 

194 

•7030 

-24.9A33 

12.83 

I9b 

6000 

198 

.7000 

-3*1 OOO 

1.56 

2ia 

6000 

223 

.9000 

-5.3000 

2.37 

?7d 

»000 

281 

.0000 

-8.2000 

2.92 

3D‘* 

6999 

306 

.1001 

- . 2002 

•07 

294 

2400 

30^' 

.3959 

-14.1599 

4*59 

3oy 

39=»9 

306 

•6031 

2.7398 

.91 

30a 

1001 

317 

.1001 

-14.0000 

4.4? 

29/ 

6999 

303 

.1001 

-5*2002 

1.72 

3oa 

3999 

3^7 

.2000 

1.2000 

.39 

29b 

bOOO 

308 

.6001 

-12.1301 

3.92 

294 

0000 

293 

.3030 

.7000 

.24 

299 

2000 

301 

.3999 

-2.2000 

.73 

3oy 

0000 

306 

.8999 

2.1001 

• 68 

31U 

B999 

309 

.3999 

1.5000 

• 48 

301 

2000 

3 '13 

.3999 

-2.2000 

. .73 

296 

2000 

291 

.3999 

4 . 8000 

1*65 

272 

/ooo 

267 

.6031 

5.0999 

1.91 

25b 

3000 

244 

.b03C 

10*8000 

4.42 

229 

9800 

231 

.0000 

-1.0200 

.44 

AVG. 

OF  HOURLY 

percent  ERRBPSa  2.79 

percent  EPR9R 

FWR 

TOTAL 

Daily  power*  101.63 
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lESOyAIQBLElQPY 


PRtDlCT  If'NI  F3R 


TUt  9 mar  76 


*««» 


PRtDlCTieN  USES  FKI 


5 Mar  76 


AS  Base  3ay« 


gain  based  9N  T-IIJK  4 MAR  76 


GAIN  IS  fixed  less  THAN  I.5, 


PRtDlCTlPN  actual  MW D I FFERENCE 


21U»!^>000 
?0l *0000 

231.4000 

221 .4000 

-20*9000 

-20*4000 

9*03 

9.21 

20S»6000 

216.9000 

-13*3000 

6.13 

194.70CO 

220*2000 

-25.5000 

11*58 

19A.U000 

224.4000 

-3Q.4000 

13.55 

?1 J»3000 

224.9000 

-11*4000 

5.07 

248.2000 

281 .8000 

-33*6001 

11*92 

298 .1001 

316.0000 

-22*8999 

7.25 

30t>»8999 

338.8999 

-32*0000 

9.46 

296.2000 

326. iQOl 

-28.9001 

8.89 

306.6997 

320.1001 

-13.4004 

4,19 

297.4331 

316.8999 

-18.9668 

5.99 

291 .2708 

294,8030 

-3.5293 

1*20 

300*2478 

296.8000 

3.9478 

1.33 

288.2998 

278.e>000 

4.7998 

1*72 

261*1543 

284 .6001 

-23.4458 

8.24 

284. 1799 

287,6001 

-3.4202 

L<l19 

291 *2000 
299.5259 
306*5000 

294 . 2000 
312.8999 
303,0030 

-3*0000 

-12*8740 

3.5000 

1*02 

4*12 

1.16 

310*7329 

288 . 7000 

22*0330 

7.63 

272.3130 

288.7000 

-11.3870 

'^•Ol 

258.9531 

258.8999 

.5532 

• 21 

252*5692 

226.7000 

25.8692 

11*41 

AVQ»  oF  HpURLV  percent  ERRf^RS*  6*06 

percent  ErRWR  F9R  T3TAL  DAIlY  P9WER»  104.21 


BESLAVAIUBLLCOBf 

fOt 


PRLJICTieM  F9R  TUt  16  “lAR  76  **** 

PRLDICTIRN  USES  FKI  12  MAR  76  AS  BASF  DAY* 


GAIN  3ASE:D  9N  thur  ^ mar  76 


GAIN  IS  fixed 

LESS  than 

1 • 5 • 

prldictipn 

ACTUAL  MW 

OIFFeRENCi 

211 . 7000 

20y»90DC 

2*1000 

20^3*2000 

190.^000 

IA.8000 

20^*^000 

19'*.1000 

11*2000 

l99. 77«9 

193.1030 

2.6789 

Ell'^+OOO 

l9S.b000 

22U.U000 

211  .«>000 

8.4000 

3.97 

2^2.3200 

2A2.6000 

-*2800 

*12 

296.6477 

302*0000 

-5.3523 

1*77 

301 *1360 

■ 3o9.booo 

-8 . 365o 

2*7o 

299.0000 

319.3000 

-20*3000 

6 . 36 

283.9000 

31 o.booo 

-27.0000 

8*70 

28/ . /QOO 

309.8999 

-22*2000 

7.16 

282*8000 

31 9*  7000 

-36.8999 

11.54 

300*3501 

306*0000 

-5.6499 

1.85 

274  * 3000 

?9b. 7o00 

-21*3999 

7.24 

279.7871 

303.6001 

-25.8130 

8.45 

277.8572 

308.8999 

-31*0427 

10.05 

296.6165 

311*0000 

-14. 383b 

4 

301*2000 

313.6001 

-l2.4ool 

3 

306*3999 

314.3000 

-7»9001 

2 

301*3401 

297.8999 

3.4402 

289.3000 

275.3000 

14.0000 

259.9571 

253. loOC 

6.4571 

244.4000 

240.0000 

4.4000 

avG.  of  hrurlv 

percent  ERRPRS= 

4.87 

PERCE-NT  ErRBR  F9R  TOTAL  DAILY  p9wER«  102. A1 


BBTjyillABliXQBf 
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ICT  1 F-R  ^aT  80  76  **«* 

PPbCMCll»M  U-cf:  sat  13  ■^^aP  76  AS  OaY. 

GAIN  Based  b\  sjn  i**  ^ia^  7t 

GA]!'.  IS  PiXtD  LE^-B  Than  l.S* 


D 1 CT  I ^•■■1  Ac1i;A|_  ^-W pI  rPEPF-^c.r 


li?/  *(’000 

iXTi.fcpou 

1 9S  t A-r) jr. 

1 6 A . 7 0 0 r 

^ • Pouo 
i5*y:oo 

?. 

9. 

1 7-d  .TQO'J 

’”^*PDj. 

1 D*  -^000 

b • 

1 . 1G74 

1 • ^DD.i 

-1  ‘BtEb 

m 

135«o=S7i 

. if  • .j  1.  y 9 

?•( 

]o/  •‘+000 

^ *3000 

5 

1 i/ • 7oqo 

1 7 A . I 0 D 0 

9»fcDJC) 

5 

PD7.U000 

^"'7.  sod: 

-•E.OOO 

EPb 

•7000 

PE-  . 7o_,o 

•0000 

? j<tj 

•^000 

■ ^P7.1oj., 

1 1 • ‘+000 

5 

c 46 

b 000 

• '^DO 

8 1 *7000 

a 

• dOGu 

EEl  •700.. 

ED' ^COO 

q 

?3i  • --*000 

= 1?^6y^0. 

1?*7000 

X 

EEE*00C0 

l^^'^DOO 

q 

21 1 • 0000 

1 

7 1 *7000 

11  • ?- 

2 10 •2000 

1 -‘.‘+00, 

23  • >■  'J0(, 

I?.b7 

21 ^ *1000 

E"  •Eooo 

15.u3 

Hi 

1 7 • 6 1 

9.95 

? • 1 r 

31  •’'  7DDo 

IP*  3000 

fi  • 7 A 

1 4 2 • y Q j r 

12*3D00 

b*lD 

lb]  • b 0 0 J 

1 R • 2b!)!) 

10*,'P 

Iho^^ono 

1 bD.SQD- 

LD'OQOO 

1 2 • *4  4 

AV"i^  oc-  ho'JPLV 

PEACE \T  tPAb^s 

= 7^7t. 

PEPOL.OT 

cA.;  Tbl  AL  D''  ! L7 

P?)-.Lb=  9j-1'j 

Ji  ‘a;  U^  r;  fu;  U/ rT  w iii 
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BKT  AVAILABLE  COPY 


PHtOlCTIP'N  F?P  SAT  ?7  '"AP  76 

PPLL1ICTIe^  USES  SAT  20  -^lAg  7fe  AS  dASE  0/Y. 
GAIN  based  0\'  SJN  cl  76 


bAi.N  IS  LESS  l-'AN  1.5* 

PRtDlCTI^N  actual  DI'FEPE^^CE 


I9t>.  AOOO 

l6A.4noc 

32*0000 

19*  A6 

17E. 700U 

158.30:3 

1 A. 2000 

8*95 

1"5'3 » ^UlTU 

I4».4o0j 

13. 8000 

13.  34 

16/ • /?CU 

1 *♦  7 . 3 0 3 0 

19. OOOO 

13.46 

lbB.;:^OCU 

1 5 1 . ‘♦OOC 

lA.boOO 

? . 5« 

lo  1 • 

1 bt . /CDD 

n . c’oou 

3TT2 

181 . IGOO 

173. 7000 

7.  AooO 

4.25 

PIU‘^000 

I'^b.lODC 

1 A. 7300 

7*50 

?L’8  .ycou  ^ 

ToFT^oco 

22*5300 

lO*9o 

?3u.lOOU 

2 1 ? . b 0 n 

17.6000 

8.2? 

?E7 .8000 

Pi::'. 603:; 

I'*  .2  000 

fi  • b5 

?'d‘*»'fUOU 

20/ • /QOO 

17.2000 

8*28 

211 *6000 

?C)0*!300C 

11*1 000 

5*54 

20i^*'*000 

1 S’:  •53000 

11*^000 

5.25 

lAy.EAoC) 

1 80*5:'00U 

b * 3^00 

4.61 

18y.l6uc 

1 7b . 7o0C 

13.A600 

7.65 

183.6900 

1 83.8033 

-•1100 

. 06 

lyy. 1300 

TFTTTooc 

1 *3?00 

.57 

Piy*uoco 

211.5033 

7*5000 

3.55 

222.1500 

21 6.800c 

5*3500 

2.47 

A . b 9 3 7 

2.20 

R*o2U0 

4.52 

8 . ES'OO 

4.78 

ltu.e857 

157.0000 

S.'S??? 

2.35 

AVG.  OF  HOURLY 

‘^ERCENT  eRR^^RS 

= 6*61 

percent  EPRBK  PBR  TOTAL  DAILY  POwfcRx  93.2? 


PRtDlCTl-^N  7t  **** 

PWLDICTI^Nl  iJSES  SUN  ?!  ^AR  76  AS  BARE 
GAIN  Based  b\  sat  ho  '^^'R  7b 


GAIK  IS  t^IXEO  LESS  T-AN  1,5. 


PRbOICTI«\' 

ACTUAL  ’'W 

01  c-ference 

W 

/• 

14u»»C00 

13B.6000 

1*2000 

*36 

136.1000 

131.3000 

4 . 6000 

3*66 

121*3000 

1 3'J*  3 000 

-8.6000 

6*61 

1 2U*3000 

1 23.7000 

-3.4000 

2.75 

121 •fiOOO 

127.0000 

-5.2000 

4*09 

IP**  *2000 

1 ?S • Sooo 

-4 .6000 

3.57 

13/.2000 

1 .tfoOD 

-7.6000 

5*25 

IbA .9ooq 

iSr.booo 

-.6000 

.39 

166  • ‘♦OOU 

177.7000 

-9*3000 

5.2? 

176.1000 

1 ?S.boot 

-13.4000 

7*d7 

1 72 . Sooo 

1 9“ .2000 

-25* A0OO 

1 ? * 5 2 

10*20 

*24 

3.31 

1 ay . 2P00 

19] .to 000 

-?.3b00 

1*23 

lS/.y600 

1 a'J.uOOC 

7.9600 

4*4? 

cO^* J400 

iSB.booO 

l5*4itoO 

8*14 

20^*Ui500 

197.6000 

7*4800 

3.79 

21i^*b000 

HC'I'  .7000 

9*8000 

4*70 

231 •boOO 

P19.2000 

1?*3qoO 

5*61 

?! j .aa  /b 

220. '♦000 

-6*5] ?b 

?*95 

20^*0200 

?C-2.boOG 

6.bc00 

3*2? 

1 86 . 6S0O 

1E9,8ooC 

-2*9200 

1 * 54 

161 *6000 

169.HO00 

-8*?0OO 

4*83 

AVG.  ttF  HPURL>^  percent  ERPBRS=  A. 44 
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vsjmMMi 


PRtUICTI9\  'iftN'  5 APR  76 

PRLDICTI^':  J5ES  P9  -^aP  76  AS  6A6E  Oay. 


GAIN  based  ev  FPI  P APR  76 
'GA  IT]— rs-p  rxT'D  L L b S 1 -i'a  M ' ] 5"* 


pPtDJCTlBN  ACTUAL 


DIPFEPE^iCE 


iBB. 1000 
iSeJ^OCOO 


0( 

15: 


1 56  *5000 
1 57 . 700C 


C'J  fT. 
Xi  to 


Til 


PKlDICTIpM 

predict  I 

F?R  IhUR  S , 

USES  inED  7 

APR  76  *•«* 

APR  76  AS  Base  D/v, 

GAIN  Based 

9M 

TJE  b APR  76 

GAIN  IS  rlxEO 

LESS  than 

pRtDlCTievi 

ACTUAL  :‘W 

X 

173. Boon 

isri.ycoo 

-12‘600U 

6.78 

I6fc.62»6 

1 7b.bC00 

* 6 * El  7 1 A 

5*?2 

IbU'BbAb 

l7l .3000 

- 1 1 » oA  bB 

6 . A 5 

lb«. 1000 

Ibb.aoOG 

2*9000 

1.39 

IbB./fce? 

lE  0*^00  0 

“11*0333 

6.10 

1 /b.200C 

190*700C 

“1  A .5000 

7.60 

213.2333 

21^+  .9000 

- 1 • 6667 

.78 

?52.277i5 

270.bo01 

“1?.3223 

6.77 

TTi  . 'o^^y 

27*'. 

-7*  *'000 

? ' 

c63.bAp7 

262.C00C' 

1 » 5 A27 

?6u • 0667 

PfrP*  7odO 

-9.6333 

3 , 

?5a  .6000 

256.0000 

-1 *2000 

.47 

?Ab.bb36 

2A3.3000 

2.5536 

1 *05 

246 . /75o 

HA9. Aoon 

-2.6250 

1*05 

2 Ay. A 000 

238.9000 

1 *5C00 

.63 

235.7667 

23A.4000 

1 .3667 

• 5H 

235.6000 

25o*Aooo 

-1  A. 6 000 

5*8? 

2A6.2000 

250*1000 

-6*9000 

2.76 

2Ab.30OO 

257.8000 

-l2*50Ul 

A. 55 

265.7312 

263.1001 

2*fclll 

• 99 

269.3000 

266*2000 

3*1001 

1*16 

235  .bA28 

PAA  , 4qq,- 

-8.8571 

3*62 

229*5857 

?2‘'  .bOOC 

.6857 

.39 

20B*b000 

l97 . 7000 

A. 9000 

?7a8 

AVG.  BF  HPURLT 

percent  erpqps 

s 3*03 

P^RLfc.'^T  FSR  total  Daily  powers  102.20 


CjCJ 

I <11  lO  lO 


BESLAyAllABLEJOPY 


PRL3ICTI9N  F3R 


ThU.^  i3  may  76 


PRLU1CTI9N  USES  12  may  76  AS  BASE  3Ay« 


gain  based  on 

T'JE  11  'IAY 

76 

GAIN  IS  FIXED 

LESS  than  1. 

5. 

PREDICTION 

actual 

difference 

% 

l<f9»AOOO 

1 59.1000 

-9‘ 7000 

6.10 

13y»9000 

145.6000 

-5.7000 

3.91 

l3a»ooou 

1 45.8000 

-7. 8000 

5*35 

l<tb.400O 

144.8000 

.6000 

.41 

13b. 1000 

147.1 000 

-12*0C0C 

8.16 

147.8800 

162.5000 

-14.6200 

9.00 

173.2000 

184.9000 

-11.7000 

6.33 

P36.b375 

2.^?*-«COO 

-6.2625 

2*58 

?5a.C356 

242.2000 

15.8356 

6.54 

24b.A400 

241 .1000 

4.3400 

1.80 

' ?6U*y7«0 

253.1000 

7.5780 

3.11 

rbl.UOOO 

237.6000 

13.^^000 

5*64 

?3U.»913 

223.9000 

6.9913 

3.12 

?34.b3b3 

,222_*b000,  . 

13.7353 

6*22 

226.9800 

222 . 8000 

4 . 1 800 

1.38 

222.1500 

209 . 2000 

12.9500 

6.19 

228.7200 

21  .0000 

14.7200 

6.88 

230*6000 

210.7000 

19.PCO0 

■ 9.44 

218.A200 

202.9000 

15.5200 

7*65 

210*  if  625 

190.1000 

20*1625 

10*61 

219.7059 

213.5000 

. 6.?05S 

2.91 

21B*7400 

220.7000 

-1.9600 

.89 

189.1059 

197,6000 

-8.4941 

4.30 

; 163.2727 

1 

169.7000 

- 6 * ** 2 73 

3.79 

1 AVG.  OF  HOURLY 

PERCENT  ERRORS*  5.12 

I 


percent  ERR3R  F9R  TOTAL  DAILY  P9WE5» 


9fa.53 


BEST jVAIUBlf. COPY 

115 

PRLDICTIPN  F3R  FRI  ?8  >'AY  76  «*** 

FRLPICTIPN  USES  THUR  27  ^AV  76  AS  BaG£ 


GAIN  Based  3N  '/<ED  =6  may  76 


gain  is  FIXE3 

less  Than 

1.5. 

PRt.DlCTI9N 

actual  MW 

Ol-FERENcE 

X 

15U»B600 

142.2000 

B.C600 

5.67 

151 *3750 

140.100C 

11.4750 

8.19 

13‘flOOO 

129.4000 

4.7000 

3.63 

133.8000 

123.3000 

10*5000 

8.52 

1 .0000 

134.3000 

9.7000 

7.22 

14y  .3000 

140.0000 

9.5000 

6.79 

I7b.y000 

168.7000 

7.2000 

4.27 

?33.?300 

222.8000 

10.4500 

4.5s 

?A8.y000 

236.0000 

12*9000 

5.47 

?A1  •‘*000 

231 .8000 

9. hood 

4*14 

24 J. 4 750 

237.3000 

6.1750 

2.60 

241.3727 

230.8000 

10*5727 

4.53 

217.3818 

227.1000 

-9.5182 

4.19 

Pll.lOOO 

221.1000 

-10.0000 

4.52 

28y.U077 

220*7000 

S.3077 

3.76 

2oi>»yooo 

217.OO0C 

-11.1 000 

5.12 

228.3233 

228.2000 

.3235 

.14 

20‘T  *829^ 

226.1000 

-21.2706 

9.41 

211.4706 

203.1000 

8.3706 

4.12 

193.3438 

200.^000 

-7.0562 

3.5? 

198.3000 

211.2000 

-12.9000 

6.11 

217.3623 

214.1000 

3.2625 

1.52 

189.7500 

188.2000 

1.3500 

.82 

154.2000 

168.4000 

-14.2000 

8.43 

AVQ.  OF  HOURLY 

percent  ERRORS'  4*89 

PERCENT  ERROR 

F8R  Tt?TAL 

daily  P9KER*  99.01 
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PRtDlCTlQN  F3R  *^EU  ?3  JUN 

PRLDICTIPN  U5FS  TUE  22  JUN 

7^  ♦**♦ 

76  AS  BASF  Day. 

GAIN  based  On 

M3N  21  JUN  76 

GAIN  IS  FIXED 

LESS  THAN  1,5. 

prediction 

actual  MW 

DICFERENtF 

X 

isy.iooo 

160.5000 

-1.4000 

.87 

147.A000 

150.2000 

-2.8000 

1.S6 

148.1000 

142.2000 

5.9000 

4.15 

141 .7000 

140.9000 

.8000 

.57 

137.8000 

136.5000 

1.3000 

.95 

147.b000 

137.1000 

10.^000 

7.59 

16y.l000 

168.4000 

.7000 

.42 

22y. 7000 

211 .0000 

1 S. 7000 

8.86 

254.6000 

249.7000 

4*9000 

1.96 

253.5600 

248.5000 

5*0600 

2.04 

266.3599 

258.3000 

8*0598 

3.12 

266.4697 

260.3000 

6.1697 

2.37 

246.0000 

251 .7000 

-5*7  000 

2.26 

253.8997 

264.8000 

-lo-yoo^ 

4.1? 

256.8000 

255.3000 

1*6000 

.59 

?4U»2000 

2*^7.0000 

-6*8000 

2.75 

247.0500 

243.1000 

3.9500 

1.62 

246.9100 

250.5000 

-3.5900 

1.43 

222.7900 

230*6000 

-7*8100 

3.39 

225.5700 

223.8000 

1*7700 

.79 

21!3»3oOO 

218.7000 

-3*4000 

1.55 

226.3500 

243.1000 

-16.7500 

6.89 

199.3000 

215.1000 

-15.8000 

7.35 

175.6000 

193.6000 

-18*0000 

9.30 

AVG*  OF  HOURLY 

percent  ERRORS 

« 3.20 

percent  EpROR  FSR  total  daily  P9WEr«  100. A7 


BESLAVAIWBlLCOa 
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PRtDICTiPN  F3R  TUE  c^  JUN  76  «*•* 

prediction  Uses  .fKl  ..  ?5  JUN  76  3asp_ 

Day*  . . 

GAIN  based  0N 

HUR  24  J'JN 

76 

GAIN  IS  FIXED 

LESS  than  1.5 

f 

PKEDICTIPN 

actual  HW 

DI.£F£ReNC£  . - 

169 .BCOO 

155.4000 

14.4000 

9.27 

153.3000 

141.9000 

11*4000 

?*.03 

144.2000 

1 27.6000 

15.6000 

15*01 

136.B00U 

121 .5000 

15.7000 

13.74 

137.3400 

122. Boon 

1 4 .4400 

-.11.75 

146./700 

129.0000 

19.7700 

15.33 

179.5000 

166.3000 

13.2000 

7.94 

20».3000 

19“. 1 000 

14.2000 

- - -.7.32 

E46. 1000 

229.6000 

16*5000 

7.19 

254.6500 

233.4000 

l9*  2500 

8.18 

260.5248 

24] . 200Q 

1 q.4248 

8.n5 

255 . 25oo 

246.3000 

8.9500 

3*63 

?51 .1000 

224.8000 

26*3000 

11.70 

P4/.6333 

.2,4  9.0000  . 

-1  * 365  7 

. *55 

23u.'34oo 

232.0000 

• 1 * 5600 

.7? 

233.8800 

232.5000 

1 *3300 

*59 

224.9462 

223 .3000 

1 .6461 

.74 

225*9600 

238.2000 

-11*2400 

4 . 72 

211 *5000 

207.1000 

4 .7000 

2.27 

20«J*5000 

2C5..300Q..  . . 

-2.5000 

1.2? 

20^*5000 

209.5000 

-6*2000 

2.96 

205*»000 

207.0000 

-1  *2000 

* 5.8 

19^*0000 

15.0.5000  - . 

lP*5nou 

6.93 

165.9000 

164.2000 

1 *7000 

1 *04 

avG*  of  hpurly 

PEHCLNT  errors*  5.14 

PEKLtNT  ErrHK  r»R  TtJT^L  DAILY  P9wER 


95.73 


beslavailable_copy 

.21 


PREDICTION 

pap  i^PI  p 

JUL 

76  ♦*** 

PREDICTION 

USES  THJA'  1 

JUL 

76  AS  base  Day. 

GAIN 

BASED 

ON  WED  3C 

JUN 

76 

GAIN 

IS  PIXLD  less  HAN 

1 • 5 

• 

PREDICTION 

actual  nw 

OIPPERENCE 

X 

1 AA 

7300 

1A7.10D0 

-2.3700 

1 

.61 

L3^ 

fapoo 

137.6000 

-A .Qoou 

2 

.91 

133 

3CU0 

1 33.b00O 

-.2D0C 

. 1 5 

131 

= 000 

130»UOOO 

1 .6000 

1 

.15 

ISA 

0333 

1H6»UOOO 

-1 .966/ 

1 

. 56 

lAE 

0000 

133.8000 

9*2000 

6 

.53 

162 

8750 

156.2000 

6.6750 

A 

•27 

?OA 

3800 

. 1S3»1000.. 

15.1P00 

8 

*03 

?3/ 

6000 

216.3000 

21*6000 

9 

.99 

?A3 

6000 

229.900Q 

13*9000 

6 

*05 

EAA 

3000 

21 7.8000 

26*5000 

12 

.17 

2AU 

1000 

227.2000 

12*9000 

5 

.68 

243 

7 000 

213.2000 

30*9000 

14 

.31 

243 

2^I>0 

22‘*.3(l00 

17.9000 

7 

.98 

230 

'♦333 

21?. 3000 

23*3333 

10 

.30 

231 

1538 

21 0.0000 

11.1538 

5 

.31 

233 

0000 

209.6000 

22*5000 

10 

.74 

233 

3867 

219.2000 

1 A. 0867 

6 

.43 

217 

9571 

206.7000 

11.P571 

5 

.45 

216 

6000 

206.3000 

10*3000 

4 

.99 

22U 

4000 

2C8.000C 

12* ACOC 

5 

.96 

22b 

6000 

197.0000 

29.6000 

1 5 

•03 

210 

0000 

187.9000 

27.1 000 

14 

.42 

179 

3000 

16?.b000 

13.8000 

8 

.34 

AVCi* 

OP  HOURLY  percent  ERR8RS=  7»09 

PFRLbNT  ERRaR  fdR  TOTAL  lJAIlY  P9WER 


93.3? 
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PKtDlCTl0\  F9P  ^EU  14  JUL  76 


PWLDICTIPN  USES  TUE  13  JUL  76  AS  BASE  JAY. 


GAIN  based  8N 

MON  12 

JUL  76 

GAIN  IS  PI^ED 

LESS  than 

1 • 5 • 

prlcictipn 

ACTUAL  MW 

DI" 

ference 

X 

13b. 6000 

135 

•3000 

.3000 

• 22 

121 .yooo 

117 

•0000 

4 

*9000 

++.19 

121*1000 

117 

*7000 

3 

*'+000 

2*89 

113*^000 

1 13 

*3000 

*6000 

• 53 

120*1000 

10J< 

•3000 

11 

*2000 

10.28 

11»*»000 

125 

*3000 

m 6 

•boOO 

5.19 

l5u*100U 

1^7 

*»000 

2 

•3o00 

1.56 

I92.b000 

19? 

*900CL_ 

. 

.+^000 

.21 

21‘**175o 

21  1 

*0000 

3 

*1750 

1 »5o 

P1»*75o6 

215 

*3000 

2 

*6500 

1*32 

217*t>000 

230 

*6000 

-13 

*?ooo 

5.72 

?22*b000 

23‘* 

*6000 

-12 

*3000 

5.24 

2iy*3ooo 

22? 

*3000 

-3 

*0000 

1.35 

226.2000 

225 

*600i3 

.6000 

.27 

250*6000 

220 

*0000 

30 

*6000 

13.91 

203*0000 

216 

*'+000 

-13 

*‘+000 

6.19 

227*7500 

214 

*9000 

12 

*8500 

5.98 

Plb.o7l4 

219 

*'+000 

-4 

*3286 

1*97 

207  *b00O 

203 

*2000 

-1 

*7000 

.81 

202  *<^000 

20** 

*1000 

-1 

*9000 

.93 

193*3000 

194 

*6000 

-1 

*3000 

.67 

Eob.aooo 

205 

*1000 

*7000 

.34 

l78*'+500 

184 

*6000 

-6 

• 3500 

3.44 

143.6000 

151 

*6000 

• s 

.0000 

5.28 

AVG.  OF  H0UWLY 

percent  ERRaRS= 

3*33 

perlent  Error 

for 

total 

daily  p9wEr=  99. 9S 

BESTAVAlUBLLCOPlf 
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PRtOlCTIPM  F9R  sat  7 AUG  76  •*** 

PRmiCTieN  USES  sat  31  JUL  76  AS  9ASF  DaY» 


GAIN  BaSEU  bn  sun  1 AUG  76 


GAIN  IS  FIXED 

LESS  T rtAN  1 1 5 . 

predict  I BN 

ACTUAL  MW 

oifferencf 

X 

I46.yooo 

151 .2000 

-4.3000 

2*84 

13U» 1000 

143.9000 

-13.8000 

9*59 

i2e;«yooo 

134.3000 

-11*4000 

8.49 

1 1000 

131.9000 

-13*8000 

10*46 

11»«0000 

126 .booo 

-8.5000 

6.72 

is^**  7600 

139,3000 

-14.5400 

10*44 

I32»y000 

141 .booo 

-B .6000 

6*08 

I4y.7000 

161.4000 

-11*7000 

7.25 

l82.b0OO 

200. 4000 

-17*9000 

8*93 

l9y .2000 

21 2 .3000 

-17.1000 

8*05 

19«»  7000 

210.9000 

-12*2000 

5*78 

197.2750 

222.4000 

-25*1250 

11*30 

lg2.«000 

211 *9000 

-29.1000 

13*73 

l8l .2857 

?06.boon 

-25.2143 

12*21 

189 .«000 

20‘**2000 

-14.4000 

7*05 

I9if  .0000 

204.100G 

-12*1000 

5.93 

193.2000 

213.9000 

-20*7000 

9.68 

SO^'^^OOO 

213.2000 

-11*8000 

5.48 

20ii»<?000 

196.6000 

5*6000 

2.85 

211 » 7ooO 

198.7000 

13*0000 

6.54 

189.0667 

H0‘**9000 

-15.8233 

7.73 

20l  't^OOC 

183.0000 

18.5000 

10.11 

191 .2000 

1 70.b000 

20*7000 

12.14 

160*4000 

152.8000 

7*6000 

4.97 

AVG.  BF  HPURLY 

percent  errprs 

» H , 10 

perc-ent  errbr 

F8R  TOTAL  Daily 

P9WER«  105.41 

PRtDICTiaN  F3R  FRl  13  AU3  76  **** 

PRLDICTIPN  USES  THUR  12  aUG  76  AS  BASE  3AY. 
GAIN  BASED  8N  W£D  11  AUG  76 
G a1  S”  fTx^  "LES^” fHAN  I.5, 


OREDICTIPN 

actual  “W 

oipference 

X 

1A4.6000 

153.1000 

-8.5000 

5*55 

I26.b667 

143.4000 

-16.8333 

11.74 

121 .6000 

137.0000 

-15.4000 

11.24 

12H .boOO 

136.1000 

-11 *6000 

8.52 

122.4000 

131.4000 

-9.0000 

6.85 

122.0000 

141 .0000 

-19.UOOO 

13.45 

lb2.O00O 

166.0000 

-lA.OOOO 

8.43 

I 9b. 4000 

210.2000 

-14.8000 

7.04 

226.8666 

240.7000 

-13.8334 

5.75 

23/. 0000 

259.8999 

-22.8999 

8.81 

24b. 9000 

256.8999 

-10.9999 

4.28 

246.4428 

263.6001 

-17.1573 

6.51 

258.3569 

260»200C 

-1.8430 

.71 

263.8000 

257.8000 

6.0000 

2.33 

26u»«999 

253.0000 

2.8999 

1.12 

254.2000 

245.1000 

9.1000 

3.71 

2A8 .3000 

253.7000 

-5.4000 

2.13 

25/. 7998 

247.2000 

10*5998 

4.29 

240.5000 

234.2000 

6.3000 

2.69 

223.2500 

232.3000 

-9.0500 

3.90 

252.2750 

238.5000 

13.7750 

5.78 

245.0000 

223.7000 

21 .3000 

9.52 

222 .1000 

200.0000 

22. 1 000 

11.05 

19U.»000 

172.5000 

18.3000 

10.61 

AVO.  of  H0URLY 

percent  errors 

* 6.50 

PERCENT  Er^OR  F0H  TOTAL  DAILY  P9WER«  101.60 
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S.QO  I.DO  12.00  19.00  10.0a  21.00 

HOUR 


B5LAVAI14B[E_C0PY 
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PRtOICTleN  F?f?  -RI  3 SEP  76  *♦** 

PPEDICTIPN)  USES  THUR  g SE°  76  A£  BASE  ^'•Y. 


GAIN  based  ^-ED  1 SEP  76 


GAIN  IS  FIXED  LESS  THAN  1,5. 

PREDlCTleN  ACTUAL  YW  DIFFERENCE 

X 

l5y»B875 

138.8000 

20.^875 

1A.76 

131. /E50 

V32.200C 

- . A750 

.36 

123»>4«33 

1 29.0000 

-5.1 167 

3. '^7 

128.1000 

1 1 8.0000 

10.1 000 

8.56 

126.8000 

132. 7000 

-6.3000 

5*20 

142.3500 

166.6000 

223.1000 

1 46.4000 

1 72.6000 
220.8000 

-4 . 0500 
-7.0000 

2.300C 

2. 

4. 

1 . 

344.6667 

237.4000 

7.S667 

3. 

261 .6635 

?44 ,6000 

1 6 . o636 

6 • 

246 .8647 

240.9030 

8*9647 

3 • 

249*2000 

236.8000 

12*4000 

5 • 

231*9556 

223.9000 

8 * 0556 

3. 

246.6000 

231 .2000 

15.?000 

6 • 

229*1000 

218. lOOC 

11.0000 

5 . 

206.6000 

215.9000 

-8*4000 

3» 

223.3667 

?12.«0D0 

10*  5667 

4 • 

223.4000 

214.0000 

9.4000 

4. 

208.6000 

1 96.9000 

11*6000 

5. 

208.6000 

206.3000 

2*6000 

1 . 

230*6000 

226.5000 

5.UOOO 

2. 

211*4000 

223.0000 

-11*6000 

5 • 

198.4615 

1 94 . 3000 

4.1615 

2. 

, 157.3286 

165.4000 

-8.0714 

4 • 

AVG*  6F  HOURLY 

PERCENT  ERRPRS 

= 4.57 

PERC-ENT  EPRHR  FbR  TOTAL  DAILY  P9wER 


97.81 


BESLAfflUSlUOa 


131 


PRhDICTlPN  FPR  SUN  5 SEP  7b  **** 


PKLDICTI9N  USES  SUN  29  aUG  76  AS  BASE  PAY* 


GAIN  Based  om 

sat  28  AUG 

76 

GAIN  IS  FIXED 

LESS  than  1 , 

5. 

PREDICTION 

ACTUAL  HW 

difference 

X 

132.9000 

128.4000 

4.5000 

3.50 

129.9000 

107.5000 

22.^000 

20*84 

1 23 . 7500 

106.6000 

17.1500 

16.09 

111.1000 

97.4000 

13.7000 

14.07 

118 .6500 

101 *0000 

17.6500 

17.48 

IIO.'ACOO 

93 . 6000 

21*8000 

23.29 

122.4000 

113.5000 

8*9000 

7.84 

130.5000. 

116.8000 

18.7000 

16*01 

14/. 1000 

144,5000 

2*6000 

1.80 

173.5000 

149.8000 

23*7000 

15.82 

177.1000 

169.2000 

7*9000 

4.67 

179.1000 

168.3000 

10*8000 

6.42 

186.6000 

169.7000 

16*9000 

9.96 

181.0000 

154.3000 

26*7noo 

17.30 

165.0000 

146.8000 

18*2000 

12.40 

169.6000 

140*8000 

28*8000 

20*45 

163.8000 

151.8000 

12*0000 

7.91 

172.8000 

164.5000 

8*3000 

5*05 

170*6000 

155.8000 

14.8000 

9.50 

174.9000 

168.5000 

6.4000 

3.80 

189.7000 

170*1000 

19*6000 

11*52 

178.9000 

161 *7000 

17.2000 

10*64 

161 *3000 

148.3000 

13.0000 

8.77 

136.4000 

127.1000 

9.3000 

7.32 

AVG*  of  H9URLY 

percent  ERRPRS*  11.35 

PERtiENT  error  for  total  DAILY  p9i^ER 


90*29 


BBCAVAIIASLE  copy 
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PREDICTION 

prediction 

F3R  PRI  10  SEP  76  ♦*** 

USFS  THJR  9 SrP  76  AS  BaSP  DaY* 

GAIN  Based 

BN 

WE 

D S SEP 

76 

gain  is  FIXED 

less  Than  i. 

r 5 • 

prediction 

actual  P'W 

^jicfe^ence 

X 

i3y«oooo 

133 

•3000 

-.9000 

•( 

128.0000 

129 

.8000 

-1.8000 

!•: 

121  •‘»000 

127 

•ROOD 

-6*1000 

4 * / 

118.1000 

126 

.0000 

-7*9000 

6*2 

121 *3000 

130 

•2000 

-8*9000 

6.5 

131 *8000 

141 

•6000 

-9.8000 

6.5 

174.4600 

177 

• «4000 

-3.3400 

1.5 

222 • 7000 

224 

•9000 

-2*2000 

psy.osob 

250 

• oooo 

-10*9500 

4.5 

P44.b000 

247 

• 9000 

-3.4000 

1.; 

P28.1000 

243 

.9000 

-17*8000 

7.5 

P36»425o 

245 

.9000 

-9.4750 

3*5 

P43. /OOO 

238 

.6000 

5.1000 

2.1 

P3J»3000 

249 

•ROOD 

-16*3000 

23/. 3000 

236 

•8000 

• boOO 

. e 

P28.yooo 

235 

•bOOO 

-6*6000 

2.5 

223. /OOO 

249 

.1000 

-25.4000 

10.: 

258.4000 

243 

•OOOO 

10*4000 

4.< 

P2b*2143 

241 

• 6000 

-15.3857 

6.: 

222*2000 

237 

•8000 

-15.1000 

6.3 

'228.4000 

239 

.5000 

-11.1000 

4.£ 

243.8000 

225 

•OOOO 

13.8000 

8.3 

18u*5000 

216 

.0000 

-35.5000 

16. A 

155*8000 

174 

•looo 

-IS. 3000 

lO.E 

Ay9*.  Qp'  HeUggLY  PE=?CE-MT  5«2g 


percent  EwRtJR  FOR  total  OAIlY  P9RLR«  104.00 


O'  O'  -r 
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PKtOlCTlBN  FSR  sat  25  SE®  76  *•** 

P'^L0ICTI9M  USES  SAT  18  SFP  76  AS  gASF  JaY. 

GAIN  Based  on  SJN  19  SEP  76 
GM N "iT^IXe'D  less  "T -IAN  1,5. 

P P lDICTION  actual  MW di^^fere'^ce % 


l6J»'t0OO 

l7j. loOO 

• 6 * 7 000 

3*94 

156.3000 

159.9000 

-3.6000 

2*25 

1 A«.5000 

iSb.BoOO 

-7.3000 

4.69 

HI  .2000 

l54,bo0O 

-13.3000 

8*61 

HI  .4000 

153.9000 

-12*5000 

8*12 

Ha*7ooo 

157 . 9ooo 

"•9*2000 

5.83 

l5a.2000 

I7o*1000 

-ll*9o00 

7*00 

l82.aoOO 

iSb.booO 

-l2*7ooO 

6*5n 

225.4000 

225.1000 

*3o00 

*13 

P3l *5000 

•231 .0000 

*5oOO 

*22 

232.5000 

231 .7000 

*9000 

*39 

233.9000 

213.3000 

20*6000 

9.66 

213.2000 

?07*5000 

5.7000 

2*75 

2H»yooo 

196.4000 

I6.bo00 

3*32 

195.9000 

191 .7000 

4.2000 

2*19 

192.4000 

183.4000 

9.0000 

4.91 

19B.7000 

193.0000 

5.7000 

2*95 

211 *1000 

199.8000 

11*3000 

5*66 

204 .»000 

190*7000 

H*1000 

7*39 

209*1750 

214*1000 

-4.9250 

2*30 

195.9000 

219.9000 

-24.0000 

10*91 

205*9000 

199.8000 

7.1000 

3*55 

171.1000 

186. 6000 

-iS.bOOO 

8*31 

147*5000 

162.2000 

-H*7000 

9*05 

AVG.  OF  HeUPLY 

percent  ERRf^RS 

» 5*23 

percent  error 

F9R  total  daily 

pSWERs  100.89 

BEST  AVAliiBLi:  COPY 


m 


PREDICTION  Ft3R  ogN  11 

prediction  uses  MON  A 

9CT  76 

5CT  76  AS  BASF  DAY. 

gain  based  on 

FRI  8 OCT  76 

GAIN  IS  FIXED 

less  than 

1.5. 

prediction 

actual  '-W 

difference 

% 

139.7000 

152.70D0 

-13.0000 

8.51 

130.7000 

143.0000 

-12.3000 

8.60 

130.3000 

135.2000 

-4.9000 

3.62 

130*3000 

137.8000 

-7.5000 

5.44 

131 .7000 

143.5000 

-11.8000 

8.22 

14/. 6000 

157 .1000 

-9.5000 

6*05 

194.7000 

201*3000 

-6.6000 

3.28 

25/.«999 

260. 6001 

-2.7002 

1.04 

277 . 7oOO 

279.6001 

-1 .9001 

• 68 

?7U»«999 

279.2000 

-S.30OO 

2.97 

264.0200 

271.1001 

-7.0801 

2.61 

264.0664 

270*7000 

-6.6335 

2.45 

256.6598 

254 .1000 

4.4598 

1.76 

260*0398 

PSS.anoo 

6.2398 

2.46 

24/. 676b 

237.6000 

10*0765 

4.24 

238.6000 

225.6000 

13*0000 

5.76 

244 . 7920 

233.0000 

11*7920 

5.06 

257.2000 

240.7000 

16.5000 

6.85 

252.9833 

246.4000 

6.5833 

2*67 

254.2429 

255.8000 

- 1 .557 1 

.61 

246.9000 

249.8000 

-2.9000 

1*16 

248.9666 

230.0000 

18.9666 

8.25 

220*2000 

218.5000 

1.7000 

.78 

198.4000 

184.8000 

13*6000 

7.36 

avg.  of  hourly 

percent  errors-  4.18 

percent  Error  ror  total  daily  power 


93.88 


BESLAVAIUBIEXOPY 
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PREDICTION  F3R  f- R I lb  3CT  76  *■>** 

prediction  U9FS  THJR  14  OCT  76  A':  BASF 

■ Day. 

GAIN  Based  on  wed  i3  9Ct  76 

GAIN  IS  FIXED  LESS  THAN  1,5. 

PREDICTIPN  ACTUAL  ^irFrRENC-  X 

• 0875 

157. 

.6286 

157. 

• /3S9- 

148. 

• 2714 

144. 

.3428 

1 55 . 

.8538 

163. 

• 2400 

193. 

*2222 

237. 

• 0498 

266 . 

•8000 

?63. 

.3398 

267 . 

.3535 

252. 

.^000 

239. 

_•  ■-'OO.Q 

24i. 

•4000 

232. 

.9818 

219. 

.9182 

231  . 

• 2030 
ii03,a_. 
•^000 
•^030 
lA-Q.D.0— 
•7000 
•^030 
1.^030.. 

• bool 
.899=3 

• 8999 


•6030 

.8030 

•lOO.G-. 

.9030 

.2030 

.J3:030 


172. *3000 


169.6030 


AvG.  of  H9U9LV  ‘^E-^CENT  E«939S  = 


0 


2.5000 

3.74 


1.71 


PREDICTION 

F4R 

PRI  12 

N9V  76  **♦* 

prediction 

USES  THJR  11 

NOV  76  AS  3ASF  DaY. 

GAIN  Based 

9N 

WED  10  N9V  76 

gain  is  FIXED 

less  Than 

1.5. 

prediction 

Actual  mw 

DlFFERENrE 

% 

isy .0000 

192.6000 

-3.1000 

1*61 

180. 1000 

190.8000 

-5.7000 

2.99 

1 70.6  000 

180.8000 

-9.2000 

4.95 

170.8000 

186.3000 

-10.5000 

5.64 

l7b.O000 

192.8000 

-16*3000 

8*45 

18/ .0000 

201«3000 

-13.8000 

6 « 86 

221.0000 

237.9000 

-16.4000 

6*89 

25/. 0000 

281.3000 

-24.3000 

8*64 

27u.U0do 

30J.2obo 

-30»2000 

10*06 

27/. 2000 

290.3999 

-13.2000 

4.55 

277.«000 

284.0000 

-6*2000 

2.18 

282.0000 

272.0000 

10*5000 

3*86 

26u.i?000 

284.2000 

-24.0000 

8*44 

267.0000 

266*5000 

.5000 

• 19 

261 .0000 

255.7000 

5.3000 

2*07 

24y .8000 

257.0000 

-7*2000 

2*80 

27**  .0243 

271.6001 

2.9248 

1*08 

302»1748 

30l *3000 

.8748 

• 29 

296.1199 

30D»2000 

-4*0801 

1*36 

292.0000 

282.3999 

10*1001 

3*58 

273.7908 

279.8000 

-6*0093 

2*15 

261.1250 

252.7000 

8.4250 

3*33 

241 .8706 

234.0000 

7.8706 

3*36 

2 2d'.  2389 

216.6000 

3.6389 

1*68 

AVG^  tiF  H90RLY  PEi^CtNT  Ef^R?RS«  

PERCENT  ERROR  FOR  TOTAL  DAILY  P9WER»  103.38 


BESLAVAILABLE  COPY 
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PREDICTION 

F9R 

THUR  1^  Ngv 

7b  «*** 

prediction 

USES  wed  17  N9V 

76  AS  BASE  DAY* 

gain  Based 

ON 

TJE  16  N9V  76 

GAIN  IS  FIXL3 

LESS  TRAN  1.5. 

prediction 

ACTUAL  MW 

difference 

X 

i9u»c:oou 

185.4000 

4*8000 

2.59 

IBI *3000 

178.3000 

3*0000 

1*68 

l7y»<i000 

176.3000 

2*9000 

1*64 

180*0000 

173.5000 

4.5000 

2.56 

iSB.yooo 

178.2000 

4 .7000 

2.64 

iy7.y5oo 

190*2000 

7 . 7500 

4.07 

23'»*8500 

227.2000 

7 . 6500 

3.37 

281 .ygBo 

277.3999 

4.5251 

1 .63 

291 *8000 

290.2000 

1*6001 

. 55 

P7O.1001 

283.3000 

-10*2000 

3.58 

?63.899y 

282.7000 

-18*8000 

6.65 

263*0000 

275.8000 

•12*8000 

4*64 

260*0000 

278.3000 

-17*8000 

6.40 

262*1001 

274.3000 

-12*2000 

4.45 

250*6000 

272.8999 

-22*2999 

8.17 

251 *8000 

269.1001 

-17*3001 

6.43 

273*8999 

290*0000 

-16*1001 

5.55 

302*0^^98 

3oft*8o0O 

-4 • 7502 

1*55 

292*3000 

296.0000 

-3*7000 

1*25 

286.b000 

290*2000 

-3.7000 

1.27 

279*3999 

280*8000 

-1*^001 

• 50 

261*0000 

256.6001 

4.3999 

1*71 

232  * 7 000 

238.8000 

-6*1000 

2.55 

20^*1571 

207*1000 

2*0571 

*99 

AVG*  tjF  HpURLY 

percent  errors 

» 3*18 

percent  Error  for  total  daily  power*  ioi»69 


b::slavailadl[  copy 
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PREDICT  I ?N 

24 

P -RI  24  OFC 

PRtD  ICTI®^ 

U5 

2S  TH'JR  23  DrC 

GAIN 

Based 

ON 

'^-D  c?  DEC  76 

GTSTTn 

IS  rlxEo 

LE 1*5* 

PRLD 

ICTI’i'N 

'’W 

difference 

% 

S'*!* 

^000 

?l'»*bOOO 

26.7(300 

12*45 

23U* 

0000 

l97.b00C 

32*4000 

16*40 

21?^ 

uooo 

^000 

S'*  * 3ooc 

12*48 

?!**  • 

0000 

1 93.4030 

1 E*6oOO 

9.5? 

21fa. 

Bcoc 

1 5'J.  7000 

EB.6300 

13*4? 

21^* 

1000 

197.0000 

22*1 000 

11*22 

24y. 

»000 

31  7.0000 

32*8000 

15.1? 

?98* 

i^ooo 

P3b.booo 

-5*7ooO 

23.55 

p^TT 

= 000 

263*2030 

32*3300 

12*15 

31b . 

OOOG 

263* iooo 

50*  7000 

1?*11 

303* 

ttOOU 

?o6  . S 03 

37*0000 

13. S7 

Tooo 

2b6 .6001 

45.5995 

17.77 

c93* 

359? 

?S3. 6003 

37.7999 

14.7° 

27/. 

tool 

252  « 6o0O 

25*0001 

9.90 

27  b"* 

ucco 

226 . ?oo  3 

49*8000 

22. 

274  . 

= 000 

24?. 4030 

31*1 000 

12.78 

?yu* 

8955 

20**  .boou 

26*3995 

9*98 

52tF 

To5o 

296.3999 

?4*3oyi 

B*2C 

3ot5* 

3995 

275.3999 

30*0000 

10*59 

310* 

^000 

2&3.0000 

30*300G 

10.82 

"7^00 

?b3.ooOC 

44.7300 

17.53 

281  * 

icoi 

?6/ * lOOl 

1 ‘‘*0000 

5.24 

261* 

ocuo 

?’b3.yo30 

10*1  cue 

4.03 

TPr, 

■^000 

231.0000 

3*6000 

1.56 

avg* 

OP  HPURLY 

percent  errors 

* 12*70 

PERt-tNT  ERPyR  TtJTAL  JAIlY  pHwtP*  88.80 


,Sl7AVAHA3lE  COPY 


• « * * 


147 


PRE.C)ICT1?N  F?R  23  DEC  76 

PWLDlCTIPN  Ur^ES  TUE 2 8_il£C _74 


AS  jAbE  day. 


gain  based  0N  mSM  ?7  DEC  76 


GAIN  IS  fixed 

LESS  Than  i # 5 . 

PRtDICTle*N 

ACTUAL  -w 

difff^rence 

% 

25t)»8oOO 

25l . 3000 

4.5000 

1.79 

2^»*6001 

242.8000 

5. 800I 

2.39 

23b  *8000 

233.2000 

2. 60OU 

1*11 

248.8000 

2^2/. 2000 

16.6001 

7.31 

254 .aooO 

238.9000 

2n.9nol 

8.94 

?5b.8999 

240.8000 

15.5999 

6.49 

280*8999 

271.8000 

9.0999 

3.35 

320*1001 

327 .7000 

-7.5999 

2.32 

35b. 8000 

341,0000 

14.8000 

4.19 

357.8999 

343,8999 

12*0000 

3.47 

368.8000 

360,1001 

3.7noo 

1 . ->9 

35o*8o00 

359,2000 

-8.8999 

2*48 

341.8000 

356.8999 

-15.5999 

4.37 

346.8000 

345.2000 

1.6001 

*Jt6 

340*8999 

338.1001 

2.7998 

.63 

334 .8000 

338.8999 

•9001 

.27 

36*8.8000 

366.8000 

-4 . 0000 

1.09 

380*8000 

368.0000 

Z> 

0 

0 

m 

r: 

4.77 

37b. lool 

365.0000 

10*1001 

2.77 

354.b403 

355.0000 

-.4597 

.13 

352.6997 

344.2000 

?.4998 

2.47 

327.0000 

341 .2000 

-14.2000 

4.16 

304. 166b 

299.0000 

5.1665 

1.73 

271.2612 

288.8000 

-17.5388 

6 . o7 

AVe*  »F  HOURLY 

PERCLnT  ERRORS 

= 3. 08 

percent  Error 

F0R  total  daily 

P9WER»  98.91 

i 
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APPENDIX  n 


FORECASTING  COMPARISONS 


Comparison  forecasts  are  made  with  the  hybird  model  and  the 
constrained  gain  method  for  selected  days. 


Affects  of  temperature  forecast  errors  on  the  load  forecast 


are  shown. 


450 — 


PHLD  1 CT  I '’N 

pop  THUP  EO 

Jan  77  *.♦* 

RKhDlCTIPN 

JJP. 

e:5 l_q 

LAN  7 7 AP 

Day, 

i 

J 

L 

GAIN  Based 

TJE  1»  JAN  77 

1 

gain  is  fixed 

LESS  than 

1 • =5  • 

1 

i_ 

PPI-.D1  CT  I 

ACTUAL  ''N 

-)  [ rpp  Pt 

y. 

1 

i 

?53.(.)1E6 

233.yODC 

l9.1l2b 

5*17 

1 

?3fa .96P6 

235.3000 

1 1 • fc  6 

B . 1 P, 

ti 

1 

E31 • /5oo 

319.10D:j 

1 2 . bBoO 

B.77 

1 

?A1  .h4itb 

2 1 b . 3 Q 0 M 

0 5 . b A 4 b 

11*51 

r?b  .O^r  l 

P)  “ i 40-',-. 

R . 1 0 Q 1 

2 . 7^ 

■f” 

23»»0000 

331.7000 

6.5000 

2.72 

?8b*(J000 

273.0000 

1 2»  OOOU 

4.40 

33/ • /onn 

31  4 *800^. 

EP.8999 

7_._2  7 

34b. b-, 00 

333.7000 

1 2.8000 

3*34 

332 • /OOO 

325.0001 

7.U999 

2.13 

__32iL?_/jD0II 

322.399- 

A • 3 0 0 0 

1 . = P 

31b. iDoi 

313. /oo  ' 

P.Ayol 

. 77 

3o« ’Pood 

3ol .399° 

7.1001 

2.36 

3 ? 3 . b 0 0 1 

3 o 2 . 7 0 0 • 1 

? T • 9 n 1 

6.-9 — 

?91 . /oOO 

297.0000 

■c  xf  " vV  * 

-5.3000 

1 .7» 

30J*3ddo 

233 . bOOO 

1 9.b0OU 

b.B7 

313  •8.-\nn 

239.89=° 

P3.9onl 

_S-»2'4 

33/ *2000 

32b. 39?^ 

1 0»80DO 

3*31 

332.bQol 

332.OO0C 

• bool 

.!« 

331  'b;?!;!! 

323 . boon 

7 . 8 D^n 

2.41 

313.8999 

312.2000 

1 *7000 

30i  ‘i^COO 

3o^»PD0rj 

-1 .3000 

?7«.Zrio<) 

Pfab  . bO"'  1 

2 p . . 1 p c?  H 

4 A-O  4 

35/ •b4b3 

?A1 .0000 

1 b . b 4 53 

6.91 

0"  riPJRL'L-PE-tCt.NX-£S--i»^PS« 


pekllnt  e^-r-sh  f?r  Total  daily  powef* 


9a  • 3D 


PKtL)lCTI'’N  F'-lR  THUR  ?0  Jam  77 


GAIN  Based  qn  tje  ii<  jan  77 


GAiNb  CU‘’PJTED  BY  SPSS  FHP  JH  ?F  PReDICTIJM 


LRt.DJ  CXJ-2-N AC.UXAt_-  ^ 


?3/  •i^OOO 
/S7-a- 
?Sy*B7oo 
?AJ.y8oo 


?4&.yA2o 
?8i ‘ysss 

_3b4 

Sbc!  • y^oo 

3AS  .6199 

31^5-5-41338- 

314 . /Aqq 

30/  •'B?0'J 

— 3gb-.b.SA^ 
296 .6879 

?0B  • /400 

3 1 3 « J 1 9-8- 

33/  .')398 
33J .8799 

?96 .»399 

313.6599 
30/. 6318 


/400 


P33.yooc, 

— 

Pl^'lOO'j 

pi^>*3o:o 



? 3 1 • 7 0 j 0 
373. 000'^ 

3 1 ‘*-  ^5^34^ 

?33.7ooj 

3?B . 6 J3 1 

IPJ- .-i39_l 

313.700,:, 
3ol  .39?=> 

7.00..0 — 

P9  7.0001; 
383 . boOO 

4?89-.  839.3 

336.3993 

333.0000 


3U’*200C 

DoP.booc 

-?bb-»6^U04 

341.0000 


-O-j-FF-r-RE-'iCE 


DOOO 
-<lEZO— - 
B700 
3*^00 

3430 

9399 

.32-98 

3600 

oiy« 

-y-ool — 

0400 

8301 


4131 

9399 

44Sli_ 

6399 

8799 

83a9- 

4600 

1318 

1396. 

7400 


1.54 

4 .69 
13.47 
— A «-2.s__ 

5. 55 

3.37 

-Irl- 54- 

5.77 

6.15 
- W4- 

.33 
1 . 93 

7-.-S6 

. 1 A 
6.6"^ 
_S.-77_. 

3.35 

.57 

— -3,.5-9- 

.47 

1 *70 

4 . a6 


;F  Haii-Ri.^ 


PE^^LtNT  EPRhR  F=)R  TrtTAL  OAIl.y  pf^JER® 


HOUR 

Constrained  Gain  Prediction  24  Jan. 


1977 
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PRLDICTIf'J  F^R  '^3N  Ja\)  77  **** 


?Rl-DirTT«M 

.ipc-p yHN 17 

1 AM_7  7 AC  =(AC--  -\AV. 

i 

GAIN  based 

1 

3N 

FRI 

21  JAN  77 

1 GAIN  IS  fixed 

i 

LESS 

Than 

1 . q . 

pRf- :)i  ri  T '»N 

APT  ! )A!  --’W 

•>  I rrrCc-Mfr 

^^*^.6333 

9000 

15.7333 

7.53 

i . 700p 

?..n‘* . 

?ooo 

^ . S c 0 0 

(1 . 1 f 

20*^  *5575 

1 96  • 

2000 

S.387b 

4. =7 

?CR»«^5oO' 

196. 

^OJC 

5.7500 

4 * <+5 

Pno . n'-i^^o 

-4-9''. 

4 '1  -1 

R . <>  'inn 

Zf  . 0 Q 

211. 

lOOG 

12.3000 

5.53 

273 . 7oOO 

262. 

1001 

11.5999 

4.43 

33  *>-.-1^1- 

_31-Z^ 

16*4  nn 1 

^ 1 ^ 

343.^300 

33  6. 

-^000 

10*0000 

2*95 

34-+  .2000 

344, 

2000 

.0000 

•00 

P4y . tj  ^00 

P3P . 

. A n 1 

L ^00 

341 .6001 

332. 

2000 

9. 4oOl 

2.83 

1 3EU«8999 

3q8. 

SODO 

12.3999 

4 . 02 

j Qiy.Hnno 

Pr'S  . 

Snn 

1 A . qnnn 

1 30^*6001 

300. 

a99R 

3 » 70O2 

2.89 

1 313.6001 

289. 

2000 

23.3000 

3 . Ob 

323 .60'^'^ 

3n7. 

39Pq 

1 A. 1 opl 

A . 24 

35t>.6B99 

329. 

3999 

26.5000 

S • 04 

361 .6001 

317. 

lOOl 

44.5000 

14.03 

34S .3000 

PQ  6 • 

7o'd  A 

^q.AQOl 

1 ? . qq 

342 . feoOO 

299. 

1001 

43.7(000 

14.61 

32b  .(1000 

281  . 

2000 

44.8000 

15.93 

?82thi7l 

-23h. 

bQ  - 

27*0171 

1 261.3462 

23o. 

1000 

31.2462 

13.58 

AVG*  »F  HSURL’t-.-gERCLMT  ER.Roagj b^-iJXX 


PERC-EMT  ERRmR  F‘i»K  ThTAl  DAJ|_y  P^^ER 


93.57 


PRtDlCTieN  F3^  'I9N  24  JaN  77 


GAIN  3aSeD  9N  F^I  21  JAN  77 


GAINS  C9^"PUTED  5Y  SPSS  F9r  NSNTH  3F  P^e'JICTI?N 
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/iEU  7 JAW  7b  *.** 

P‘^LlUCTIPNi  tJ?F?  TU£  b JA‘-I  7d AS  tiAWP  3Ay. 

gain  SaSFD  S3\'  vt*:j  5 j^N  7A 


jjj';  is.?^ 

c2C)*A99y  P3P.P0JG  -]2*7q(j1  5»a? 

? lo  • Aif  j(j  • 2oj;i  - 1 ? t 3 Sou  6 . A 0 

P33.2A01  P?3.<^,'):0  p*aaq1  a. 54 

c'2b*AHol  ?2^.OUJ0  1,11 

?JA  .rtAni  PAif  ..-foj';-  ^mnrg^  3T^‘ 

26y,bG'J0  '3  o'j'jr  -pa.Squo  11»?? 


PE‘<'-EnT  L'?-<«K  F^P  TttiAL  JAILY  ^''^VitSs  37.53 


t58 


j prediction 

FOF 

«fO  7 JAN  76  «*#« 

: PREDICTION 

USES  TUE  6 JAN 

7b  AS  BASE  Day. 

1 

GAIN 

based 

9N 

6:3N  5 JAN  76 

GAIN 

IS  FIXED 

LESS  Than  1.5. 

phedictihn 

actual  ‘'W 

01  - FERENCE 

% 

i 2b3 

3?99 

293.300C 

-39*9001 

13*60 

■ 23b 

8999 

233.2000 

2*6999 

1*15 

232 

0000 

229.2000 

2*8000 

1*22 

23b 

8000 

223.8000 

13*0000 

5*31 

241 

1001 

228.0000 

13*1001 

5*75 

251 

6001 

242.3000 

9*3C0l 

3*84 

28U 

bOOO 

303.8000 

-23*  30OC 

7*67 

33b 

i^OOO 

33“. 3000 

*8999 

*27 

343 

3599 

337.2000 

6.20OC 

1*34 

341 

1001 

333.OO00 

8*1 OOl 

2*43 

33b 

bCOO 

337.3000 

-1*8000 

• 53 

32y 

6001 

322.8999 

6*  7002 

2*08 

313 

bl78 

303.7000 

9*8179 

3*23 

321 

2141 

30“«6001 

16*6140 

5*45 

3oy 

0000 

298.8000 

10*2000 

3*41 

30b 

2712 

281 .3000 

24*9712 

8*38 

33/ 

/66b 

312.8999 

24.8667 

7*95 

; 34y 

7856 

333.0000 

16*7856 

5*04 

‘ 33/ 

0747 

324 .1001 

12*9746 

4*00 

328 

7664 

300*8999 

27 • 8665 

9*26 

31b 

0222 

297.8999 

17*1223 

5*75 

2^/ 

6499 

278.8000 

18*8499 

6*76 

288 

1001 

248.3000 

39*8001 

16*03 

25y 

7122 

226.2000 

33*5122 

14.82 

AVCi* 

OF  HOURLY 

percent  errors 

o 

• 

u> 

II 

PEKt-ENT  ERS»R  r^R  TOTAL  DAILY  P5wER  = 


yfa.50 
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PREDlCTieN  F9R  SAT  31  JAM  76  **** 


prediction 

USfS  sat  24  JaN  76 

AS  BASF 

Day. 

GAIN  based 

BN  SJN  25  JAN 

76 

GAINS  ce^iPUTED  BY  SPSS  FOR 

NBNTH 

9F  PReDICTI?"N 

PREDICTION 

actual  pw 

djfflrencz 

X 

?7b 

1899 

245.4000 

29* 

7899 

12.14 

27u 

300C 

236.4000 

33. 

9001 

14.34 

?6b 

0898 

230»OODO 

36* 

0898 

15.69 

25U 

6600 

229.8000 

20* 

8600 

9.08 

25U 

1100 

228.9000 

21. 

2100 

9.27 

263 

5598 

231 .8000 

31  . 

7598 

13*70 

274 

5000 

257.6001 

16. 

8999 

6 . 56 

30B 

2698 

274.1001 

2S. 

1697 

/■ 

in. 28 

330 

9497 

294,5000 

36. 

4497 

/ 

12.38 

33/ 

2798 

297.6001 

39* 

6797 

/ 

13.33 

352 

7698 

298.5000 

54. 

2698 

IB.  18 

354 

2500 

3CO»500C 

53. 

7500 

17.89 

330 

8799 

297 .7000 

42. 

1799 

14.17 

337 

8799 

275.8000 

62* 

0798 

P2.51 

31» 

0641 

267.8999 

51* 

0642 

19*06 

32^ 

6240 

278.3000 

47. 

3?40 

17.00 
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284.8000 
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6399 
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10*45 
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315.5000 

32* 
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338 

3198 
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_.33. 
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in*96 

314 

0797 

300*6001 

14. 
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4.78 

290 

6001 

275.6001 

24. 
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8.71 

271 

6797 
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_ 12*. 

9797 

5.n? 
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3600 

240.7000 

8. 

6600 
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APPENDIX  5 

LOAD  FACTORS 

\ 

Central  Vermont  Public  Service  Corporation  daily  load  factors 
are  presented  for  the  first  100  days  in  both  1976  and  1977. 
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LOAD  FACTORS  FIRST  100  DAYS  1976 


DAY  LOAD  FACTOR 

MIN. 

MAX. 

DAILY  TOTAL 

1 79.80 

165.9 

291.4 

5580.59 

2 80. 5U 

204.5 

347.4 

6715.18 

3 82.43 

188.3 

305.3 

6040.09 

4 79.76 

167.5 

290.0 

5551.29 

5 81.06 

212.5 

375.5 

7305.38 

6 85.94 

254.0 

365.4 

7536.98 

7 85.59 

223.8 

337.3 

6928.78 

8 82.94 

189.8 

345.4 

6875.09 

9 86.67 

245.3 

363.0 

7551.09 

10  82.24 

234.7 

353.9 

6985.19 

11  88.29 

248.4 

322.9 

6841.78 

12  84.93 

210.3 

350.2 

7138.09 

13  82.85 

216.5 

342.2 

6803.99 

14  83.67 

174.9 

321.3 

6451.78 

15  82.87 

202.6 

345.0 

6861.68 

16  83.04 

226.7 

363.9 

7252.28 

17  79.71 

204.1 

349.5 

6686.09 

18  88.45 

236.8 

317.7 

6744.09 

19  83.31 

248.1 

388.7 

7771.48 

20  84.72 

231.3 

355.9 

7236.58 

21  85.27 

196.7 

325.5 

6663.08 

22  80.15 

191.8 

378.0 

7271.39 

23  87.58 

274.6 

386.9 

8132.38 

24  86.98 

248.4 

354.7 

7404.78 

25  85.11 

216.8 

294.9 

6023.49 

26  79.42 

190.7 

342.7 

6532.18 

27  81.47 

170.0 

3I6.6 

6190.29 

28  83.82 

180.0 

323.1 

6499.58 

29  85.21 

206.5 

324.1 

6627.69 

30  82.67 

202.0 

342.4 

6793.19 

31  86.24 

228.9 

315.5 

6530.38 

32  85.60 

186.9 

266.3 

5470.89 

33  78.27 

164.1 

361.5 

6790.69 

34  85.86 

220.8 

338.5 

6975.09 

35  86.12 

218.7 

335.1 

6925.88 

36  85.27 

228.5 

338.2 

6921.39 

37  85.35 

209.7 

333.7 

6835.68 

38  86.38 

223.2 

308.9 

6403.88 

39  87.29 

195.4 

262.3 

5494.99 

40  82.57 

181.2 

317.9 

6299.99 

41  81.78 

201.9 

321.7 

6313.89 

42  83.94 

165.5 

291.2 

5866.08 

43  85.02 

185.3 

300.9 

6139.49 

44  84.39 

172.3 

289.6 

5865.68 

45  82.41 

190.2 

305.3 

6038.38 

46  87.02 

180.1 

273.9 

5720.59 

47  79.55 

155.2 

298.4 

5697.19 

48  83.51 

168.6 

308.4 

6180.79 

49  83.66 

181.4 

317.1 

6366.79 

50  82.87 

176.6 

309.8 

6161.58 
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51  82.88 

185.2 

313.6 

6237.69 

52  87.09 

194.3 

274.7 

5741.69 

53  81.68 

149.0 

249.0 

4881.29 

54  81.62 

187.8 

342.4 

6706.89 

55  81.16 

220.7 

337.0 

6564.38 

56  82.51  ‘ 

183.8 

295.5 

5851.59 

57  83.06 

168.9 

289.4 

5769.29 

58  83.52 

165.9 

281.0 

5632.49 

59  80.99 

160.5 

283.5 

5510.39 

60  83.73 

174.3 

250.1 

5025.69 

61  80.34 

148.5 

301.6 

5815.09 

62  84.44 

180.2 

319.1 

6467.08 

63  84.80 

194.0 

325.4 

6622.19 

64  84.93 

195.3 

306.7 

6251.48 

65  83.67 

169.0 

281.4 

5650.99 

66  80.56 

160.5 

275.8 

5332.39 

67  83.55 

174.4 

255.6 

5125.39 

68  83.38 

182.5 

321.8 

6439.63 

69  81.83 

216.9 

338.4 

6645.89 

70  84.80 

201.4 

306.1 

6229.79 

71  82.14 

178. 1 

311.4 

6138.89 

72  84.53 

204.2 

312.9 

6348.18 

73  85.81 

173.9 

261.3 

5381.09 

74  85.95 

164.2 

226.0 

4661.99 

75  31.15 

153.1 

292.9 

5704.79 

76  85.07 

190.4 

319.7 

6527.18 

77  83.76 

197.4 

324.6 

6525.59 

78  86.82 

212.6 

313.4 

6530.39 

79  82.37 

193.6 

311.9 

6165.68 

80  85.83 

158.9 

227.1 

4677.89 

81  77.45 

119.3 

219.6 

4082. 10 

82  80.25 

142.0 

310.4 

5977.99 

83  81.05 

202.0 

305.4 

5940.99 

84  81.48 

165.9 

274.8 

5373.80 

85  85.15 

147.7 

259.9 

5311.09 

86  82.41 

157.6 

261.5 

5171.89 

87  85.04 

147.8 

216.8 

4424.70 

88  79.85 

123.7 

220.4 

4223.89 

89  80.07 

151.7 

280.0 

5380.79 

90  81.05 

163.0 

269.0 

5232.49 

91  81.89 

144.1 

261.7 

5143.29 

92  80.87 

146.6 

283.3 

5498.79 

93  80.75 

166.0 

279.9 

5424. 19 

94  82.89 

148.6 

237.5 

4724.69 

95  80.66 

142.3 

217.2 

4204.70 

96  78.56 

149.5 

285.1 

5375.29 

97  84.36 

153.4 

261.1 

5286.39 

98  83.56 

159.2 

266.9 

5352.48 

99  82.86 

165.8 

278.9 

5546.59 

100  82.38 

159.0 

277.1 

5478.69 

LOAD  FACTORS  FIRST  100  DAYS  1977 
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DAY  LOAD  FACTOR 

MIN. 

MAX. 

DAILY  TOTAL 

1 83.711 

217.7 

317.2 

6375.18 

2 81.96 

195.8 

288.6 

5676.98 

3 83.18 

192.1 

340.0 

6787.18 

4 83.72 

200. 1 

340.9 

6849.49 

5 84.47 

218.2 

346.8 

7030.68 

6 82.97 

227. 1 

362.6 

7219.98 

7 86.37 

202.2 

327.2 

6782.79 

3 81.54 

198.7 

323.8 

6336.99 

9 87.92 

223.6 

291.1 

6142.59 

10  80.41 

209.3 

361.1 

6969.09 

11  81.38 

206.6 

369.2 

7211.19 

12  85.10 

235.5 

365.6 

7467. 18 

13  85.81 

242.0 

363.5 

7588.69 

14  85.63 

241.0 

361.0 

7418.99 

15  85.63 

219.9 

314.4 

6461.59 

16  33.72 

197.4 

294.7 

5921.59 

17  81.84 

212.7 

385.6 

7573.38 

18  87.15 

264.4 

384.2 

8036.07 

19  84.31 

247.1 

363.5 

7355.58 

20  35.00 

216.3 

333.7 

6807. 18 

21  84.66 

208.6 

334.0 

6786.18 

22  80.59 

211.6 

344.5 

6663.59 

23  85.35 

212.8 

293.0 

6001.79 

24  80.75 

196.2 

344.2 

6670.89 

25  83.08 

189.6 

322.0 

6420.19 

26  82.81 

182.7 

331.4 

6586.48 

27  82.69 

213. 1 

351.5 

6975.69 

28  86.43 

237.8 

358.3 

7432.08 

29  82.85 

220.4 

346.2 

6883.79 

30  88.32 

232.4 

303.7 

6437.48 

31  84.26 

221.1 

365.3 

7387.68 

32  84.56 

227. 1 

359.7 

7299.49 

-3  83.63 

238.8 

364.5 

7316.39 

34  36.05 

239.4 

349.6 

7219.88 

35  84.61 

208.8 

326.8 

6635.89 

36  34.34 

204.9 

303.5 

6143.59 

37  83.14 

186.9 

292.3 

5832.19 

38  83.02 

212.6 

356.9 

7110.79 

39  84.36 

230.0 

345.7 

6998.98 

40  85.34 

215. 1 

331.0 

6779.08 

41  82.72 

189.4 

313.0 

6213.59 

42  82.90 

198.3 

315.2 

6270.88 

43  82.76 

169.3 

263.8 

5339.19 

44  83.05 

146.6 

241.6 

4815.59 

45  31.96 

150.1 

300.5 

5911.19 

46  32.32 

178.4 

319.2 

6305.99 

47  84.49 

208.4 

337.9 

6852.09 

48  86.02 

237.8 

351.4 

7254.  18 

49  83.43 

228.6 

344.8 

6903.68 

50  85.14 

202.5 

298.5 

6099. 19 
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51  81.20 

176.6 

282.2 

5499.79 

52  82.00 

178.6 

323.1 

6358.38 

53  8M.U6 

217.5 

343.3 

6953.79 

54  82.28 

198.4 

337.9 

6672.78 

55  84.81 

203.6 

332.4 

6765.68 

56  82.70 

186.7 

312.5 

6202.79 

57  84.92 

174.4 

266.5 

5431.50 

58  80.74 

159.8 

249.6 

4836.89 

59  81.95 

154.4 

301.7 

5934.19 

60  83.39 

179.1 

305.7 

6118.29 

61  84.64 

193.9 

311.5 

6327.88 

62  81.76 

189.7 

307.6 

6036.08 

63  84.08 

184.6 

312.8 

6311.99 

64  84.22 

167.1 

260.5 

5265.59 

65  82.54 

151.9 

245.3 

4859.19 

66  83.73 

157.7 

295.4 

5935.79 

67  84.34 

177.4 

287.8 

5825.49 

68  81.93 

175.3 

284.1 

5586.09 

69  84.22 

146.4 

255.9 

5172.59 

70  83.10 

153.8 

262.2 

5229.09 

71  79.68 

134.4 

238.4 

4553.69 

72  83.36 

125.9 

215.5 

4311.20 

73  80.77 

126.8 

277.4 

5377.59 

74  82.77 

155.3 

276.0 

5482.79 

75  83.05 

157.0 

286.2 

5704.69 

76  85.14 

174.7 

294.3 

6013.79 

77  85.82 

191.1 

305.8 

6298.69 

78  86.24 

182.0 

261.8 

5418.49 

79  85.02 

173.7 

236.4 

4823.89 

80  81.06 

157.3 

293.8 

5715.88 

81  81.08 

173.5 

313.6 

6102.29 

82  81.93 

159.4 

295.9 

5818.59 

83  85.03 

177.4 

297.2 

6065.29 

84  84.36 

182.4 

309.3 

6262. 19 

85  87.01 

183.5 

256.9 

5364.39 

86  87.64 

157.3 

213.7 

4494.89 

87  80.22 

142.0 

278.4 

5359.89 

88  79.97 

146.8 

272.4 

5227.98 

89  80.38 

135.3 

264.8 

5108.29 

90  80.10 

135.8 

270.9 

5207.99 

91  80.04 

160.8 

279.4 

5367.39 

92  84.11 

158.5 

244.5 

4935.49 

93  80.54 

129.5 

216.2 

4179.00 

94  81.25 

154.5 

285.4 

5564.98 

95  81.71 

160.5 

293.6 

5757.69 

96  83.76 

169.5 

284.5 

5718.89 

97  84.07 

175.7 

291.1 

5873.49 

98  83.38 

168.9 

282.2 

5647.39 

99  36.58 

180.1 

254.9 

5296.38 

100  86.99 

157.3 

211.4 

4413.70 

j 


